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1.  Introduction 

1.1.  Motivation 

Unmanned  Underwater  Vehicles  (UUV’s)  perform  many  of  their  missions  in  shallow 
water  environments  subject  to  the  forces  of  ocean  waves  and  the  proximity  to  the  ocean 
floor.  Under  these  conditions,  accurate  vertical  position  control  is  necessary  to  prevent 
broaching  or  hitting  the  ocean  floor.  Accurate  horizontal  position  control  is  necessary  to 
enable  the  UUV  to  conduct  its  mission  with  accuracy  and  return  to  a  predetermined 
recovery  point.  Shallow  water  position  control  is  made  more  difficult  by  ocean  waves. 

In  deep  water  the  effects  of  these  waves  are  negligible,  but  the  effects  in  shallow  water 
are  significant.  Important  shallow  water  missions  include  pollution  monitoring,  marine 
life  sampling,  bottom  contour  mapping,  and  mine  location. 

Currently,  UUV’s  are  controlled  in  shallow  water  by  altering  empirical  control 
parameters  for  better  shallow  water  performance  and  by  establishing  empirically  based 
operating  depth  limits  on  the  UUV  operations.  These  operating  depth  limits  are  based 
upon  wave  conditions.  With  a  thorough  understanding  of  the  dynamics  of  UUV’s  in 
shallow  water  and  the  forces  and  moments  on  vehicles  due  to  sea  waves  in  these  waters, 
improved  control  systems  and  vehicle  designs  can  be  achieved  to  allow  the  UUV  to 
operate  in  shallower  water  and  in  larger  waves  than  is  commonly  done.  This  will  allow 
the  UUV  to  be  more  effectively  perform  its  missions. 

This  thesis  explores  the  effects  of  variation  in  water  depth  and  vehicle  submergence 
on  added  mass,  damping,  and  restoring  forces. 

1.2.  Historical  Background 

This  work  builds  on  the  work  that  Timothy  Prestero  performed  to  build  a 
mathematical  simulation  of  the  REMUS  behavior  in  deep  water.  This  work  is  reported  in 
his  Master  of  Science  Dissertation  for  the  Massachusetts  Institute  of  Technology /Woods 
Hole  Oceanographic  Institution  Joint  Program  in  Oceanography/Applied  Ocean  Science 
and  Engineering  entitled  “Verification  of  a  Six-Degree  of  Freedom  Simulation  Model  for 
the  REMUS  Autonomous  Underwater  Vehicle.”  Mr.  Prestero  calculated  the 
hydrodynamic  and  hydrostatic  coefficients  based  upon  deep  water  performance  far  from 
a  boundary  surface.  This  thesis  extends  Mr.  Prestero’s  work  by  determining  those 
coefficients  in  shallow  water  and  near  the  surface. 1 

1.3.  Research  Platform 

The  platform  for  this  research  is  the  REMUS  (Remote  Environmental  Monitoring 
UnitS)  AUV  (autonomous  underwater  vehicle)  developed  by  the  Oceanographic  Systems 
Laboratory  at  the  Woods  Hole  Oceanographic  Institution.  This  low-cost,  modular  AUV 
was  developed  for  coastal  monitoring  and  multiple  vehicle  survey  operations.2  REMUS 
has  also  been  adopted  for  use  in  mine-counter  measure  operations  for  the  United  States 
Navy.3  REMUS  has  most  recently  been  used  by  the  United  States  Navy  to  hunt  for 
mines  from  the  Iraqi  port  of  Umm  Qasr  in  support  of  Operation  Enduring  Freedom.4 


1.4.  Assumptions 

To  simplify  analysis,  the  author  made  the  following  assumptions: 

•  The  vehicle  is  port-starboard  symmetric. 

•  The  vehicle  is  a  rigid  body  of  constant  mass. 

•  There  are  no  significant  vehicle  dynamics  occurring  faster  than  the  data  sampling 
frequency  of  25  Hz. 

2.  The  Coordinate  System 

The  coordinate  system  used  for  this  research  is  shown  in  Figure  1 .  This  is  a  body- 
fixed  right-handed  coordinate  system  with  the  x  axis  defined  along  the  axial  length  of  the 
vessel  and  the  z  axis  defined  downward.  The  origin  of  the  body-fixed  coordinate  system 
is  at  the  vessel  amidships.  The  variables  shown  in  Table  1  are  defined  using  the 
coordinate  system  shown  in  Figure  1 . 


Figure  1.  Sketch  showing  positive  directions  of  axes,  angles,  velocities,  forces  and  moments. 
(Feldman,  1979) 


Table  1.  Coordinate  System  Variables. 


X 

y 

z 

u 

V 

w 

p 

q 


Surge  position  forward. 

Sway  position  to  the  right. 
Heave  position  downwards. 
Velocity  in  the  surge  direction. 
Velocity  in  the  sway  direction. 
Velocity  in  the  heave  direction. 
Rotation  about  the  x  axis. 
Rotation  about  the  y  axis. 


r  Rotation  about  the  z  axis. 
X  Force  in  the  x  direction. 
Y  Force  in  the  y  direction. 

Z  Force  in  the  z  direction. 
K  Moment  about  the  x  axis. 
M  Moment  about  the  y  axis. 
N  Moment  about  the  z  axis. 


3.  The  Equations  of  Motion 

This  work  primarily  explores  the  forces  and  moments  in  sway,  heave,  pitch  and  yaw 
due  to  motion  of  the  vehicle  itself  in  finite  depth  water.  For  each  series  of  tests,  the 
vessel  was  moved  in  only  one  plane  at  a  time.  The  resulting  hydrodynamic  forces  were 
determined  by  subtracting  inertial  forces  from  the  measured  forces.  Then,  the 
hydrodynamic  coefficients  were  extracted  from  the  hydrodynamic  forces.  The  equations 
of  motion  were  used  to  perform  the  mathematical  operations. 

The  forces  and  moments  associated  with  surge  and  roll  have  not  been  investigated. 

3.1.  Vessel  Inertial  Dynamics 

The  linearized  equations  of  motion  with  a  body-fixed  coordinate  system  for  an 


unrestrained  vessel  in  water  are  given  by 

X  =  m\ 

[u-vr  +  wq- xG(q2  +  r2)  +  yG(pq  -r)  +  zG(pr  +  q)j 

Y  =  m\ 

[v-wp  +  ur-yG(r2  +  p2)  +  zG(qr-p)  +  xG(qp  +  r)j 

Z  =  m| 

\w-uq  +  vp-zG(p2  +  r2)  +  xG(rp-q)  +  yG(rp  +  p)] 

k  =  ixxp+ (4  -  -  4  (p<i  +  0 + 402  -  v2) + 4  0>r  ~^+ 

m[yG(w+pv-qu)-zG(v  +  ru- pw)]  (1) 

M  =  I„q  +  (4  -  la)pr  -  4(p  +  qr )  +  IJp2  -  r 2)  +  I^qp  -  r)  + 
m[zG(u  +  qw  —  rv)  —  xG{w  +  pv-qu )] 

N  =  Izzr  +  (Iyy-Ixx)pq-Iyz(q  +  rp)  +  Ixy(q2-p2)  +  Ixz(rq-p)  + 
m[xG(v  +  m  —  pw)  -  yG(u  +  qw  -  rv)] 

where 

m  is  the  mass  of  the  vessel 

(xo,yG,ZG)  are  the  coordinates  of  the  center  of  gravity  of  the  vessel  in  the  body 
fixed  coordinate  system. 

Ijk  are  the  moments  of  inertia. 

These  equations  can  be  simplified  by  fixing  the  coordinate  system  at  the  midship 
location  of  the  vehicle.  The  equations  can  also  be  simplified  by  assuming  that  the  lateral 
distance  from  the  midship  location  to  the  center  of  gravity  is  negligible,  i.e.  yc  =0. 

Further  simplification  can  be  obtained  by  testing  and  analyzing  motions  in  the  vertical 


and  horizontal  planes  separately.  This  research  does  not  examine  hydrodynamic  forces  in 
surge  and  roll,  so  the  relevant  simplified  equations  are: 

Y  =  m  [v  +  Ur  +  xGr\ 

Z  =  m[w-Uq  -  zGq 

M  =Iyyq  +  m[zG(u-vr  +  wq)-xG(w-Uq)] 

N  =  lj  +  mxG  (v  +  ru) 


3.2.  Hydrodynamic  and  Hydrostatic  Equations 


This  thesis  explores  the  hydrodynamic  forces  and  moments  due  to  unsteady  motion  of 
an  underwater  vehicle.  For  that  reason,  hydrostatic  effects  have  been  removed  from  the 
data  by  subtracting  the  mean  forces  and  moments  from  all  measured  forces  during  the 
analysis. 

The  forces  and  moments  experienced  by  a  ship  are  assumed  to  be  the  forces  and 
moments  arising  from  motions  of  the  ship  which  in  turn  have  been  excited  by  another 
source.  These  forces  and  moments  are  computed  as  functions  of  speed  and  acceleration. 
A  mathematically  useful  form  is  derived  using  the  Taylor  expansion  of  a  function  of 
multiple  variables.  For  example,  sway  force,  Y,  and  yaw  moment,  N,  are  represented 
functionally  as 


Y  =  Fy(u,v,u,v,r,r) 

N  =  Fr(u,v,u,v,r,r) 

The  Taylor  expansion  of  a  single  variable  states  that  if  the  function  of  a  variable,  x, 
and  all  its  derivatives  are  continuous  at  a  particular  value  xj,  then  the  value  of  the 
function  at  a  value  of  x  close  to  xj  can  be  expressed  as 


dx- 


dx 


dx2 


dim , ... ,  * 


,  dnm 


dx 3 


dx" 


(4) 


where  > 


f(x)  is  the  value  of  the  function  at  x  close  to  xj 
f(xj)  is  the  value  of  the  function  at  x  =  Xj 
bx  =  x  -  Xi 


d"f(x)  . 


dx" 


is  the  rtth  derivative  of  the  function  evaluated  at  x  =  xi 


By  making  6x  sufficiently  small,  higher  order  terms  can  be  neglected.  Equation  (4) 
reduces  to 

df(x) 


f(x)  =  f(xl)  +  Sx- 


dx 


(5) 


and  is  called  the  linearized  form  of  the  Taylor  expansion. 

For  functions  of  two  variables  the  linearized  form  of  the  Taylor  expansion  is 

ox  dy 

Again,  6x  and  6 y  must  both  be  small  enough  that  higher  order  terms  can  be  neglected. 


(6) 


Hydrostatic  motion  stability  typically  considers  the  effect  of  very  small  perturbations 
on  the  behavior  of  the  ship.  Thus,  the  linearizing  assumption  for  the  Taylor  expansion 
can  be  used  to  describe  the  hydrodynamic  behavior  of  a  body.  Analysis  of  data  from  this 
research  indicates  that  similar  non-dimensional  results  were  obtained  for  tests  done  at 
different  amplitudes  and  the  Fourier  coefficients  at  the  excitation  frequencies  dominated 
all  others.  Because  of  these  facts,  the  linear  terms  do  indeed  predominate  and  the  model 
based  on  them  is  sufficient  to  describe  the  relation  between  vehicle  motions  and  the 
forces  and  moments  they  generate.  Using  the  linearized  Taylor  expansion,  equation  (3) 
can  be  written  as 

Y  =  F(ul,vvux,vi,ri,rl)  +  (u-ul)—+(v-vl)—+...  +  (r-rl)— 

du  dv  dr 

X7  .  .  . v  /  . dN  ,  . dN  ,.  .^dN  {) 

N  =  Fr(ui,vltu^,vvr^,r^)  +  (u-u^)-—  +  (y-v^)—  +  ...  +  (r-rx)— 

du  dv  dr 

At  this  point  several  simplifying  assumptions  can  be  made.  The  first  assumption  is 
that  the  initial  motion  is  in  a  straight  line  at  some  constant  speed.  Therefore, 
n,=Vi=^=rj=0.  The  ship  is  symmetrical  about  the  xz-plane,  so  v/  =  0.  Symmetry 

also  leads  to  the  conclusion  that  dY/du  =  dY/du  =  0  because  forward  motion  will  not 
cause  a  lateral  velocity.  Also,  a  ship  traveling  forward  in  equilibrium  in  straight  line 
motion  experiences  no  sway  force,  so  the  term  Fy(ux,vx,ux,vx,rx,rx)  is  also  zero.  The  term 

uj  is  equal  to  the  straight  line  velocity  U.  These  assumptions  reduce  equation  (7)  to 

v  dY  dY  .  dY  dY  . 

dv  dv  dr  dr 

^  dN  dN  .  dN  dN  .  1  ' 

dv  dv  dr  dr 

In  the  simplified  notation  used  by  the  Society  of  Naval  Architects  and  Marine 
Engineers  and  including  Pitch  and  Heave,  the  simplified  linear  hydrodynamic  equations 
become5 

Y  =  7vv  +  y^v  +  Yrr  +  Y-r 
N  =  Nv  +  N-v  +  Nr  +  N.r 

(9) 

Z  =  Zww  +  Z^w  +  Zqq  +  Zfl 
M  -  Mww  +  M^w  +  Mqq  +  M^q 

The  simplified  notation  is  interpreted  such  that  Yvv  is  the  sway  force  related  to  sway 
motion  and  Yv  is  the  maneuvering  coefficient  of  sway  force  due  to  sway  motion. 

In  accordance  with  the  standard  notation  the  terms  of  equation  (9)  include  the  effect  of 
the  rudder  and  stem  planes  held  at  zero  degrees.  The  experiments  to  extract  the 
coefficients  were  all  performed  with  no  deflection  of  the  control  surfaces.  Other 
experiments  were  performed  with  control  surface  deflection.  For  those  experiments, 
equation  (9)  has  additional  terms  related  to  rudder  and  stem  plane  angle.6 

3.3.  Added  Mass  and  Damping 

The  hydrodynamic  forces  relating  to  the  motion  of  the  body  in  the  fluid  can  be  divided 
into  components  in  phase  with  the  acceleration  and  components  in  phase  with  the 
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velocity  of  the  body.  The  hydrodynamic  force  due  to  the  acceleration  of  the  body  in  a 
fluid  is  known  as  an  added  mass  force.  The  hydrodynamic  force  due  to  the  velocity  of 
the  body  in  the  fluid  is  known  as  a  damping  force.  These  forces  are  can  be  discerned  by 
their  phases  relative  to  the  driving  motion.  Forces  in  phase,  but  opposite  in  sign,  with  the 
driving  motion  are  related  to  acceleration  and  are  added  mass  forces.  Forces  90  degrees 
out  of  phase  with  the  driving  motion  are  related  to  velocity  and  are  damping  forces.  In 
terms  of  complex  notation,  the  added  mass  is  related  to  the  real  component  of  the 
measured  force  and  the  damping  is  related  to  the  imaginary  component  of  the  measured 
force. 

4.  Non-Dimensionalizing 

Throughout  this  thesis  several  quantities  are  given  in  both  dimensional  and  non- 
dimensional  form.  Final  results  are  given  in  non-dimensional  form  to  be  readily  available 
for  use  with  other  bodies  of  similar  shape.  Non-dimensional  quantities  are  denoted  by  a 

7 

prime  symbol  (').  The  equations  for  non-dimensionalizing  are: 

Y'  = - — - 

\pU2L2 


\pUl} 


ipuz 

Y 


\pU2L2 


\pUD 
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Other  non-dimensional  equations  include 
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Using  non-dimensional  coefficient,  the  equations  of  motion  have  the  form 

Y'  =  Y.'v'+  Y.'v'+Y  'r'+Y  'r' 


(12) 

(13) 


5.  Experimental  Procedure 

The  determination  of  the  maneuvering  coefficients  was  conducted  using  both  full 
scale  and  model  scale  experiments.  Full  scale  experiments  were  used  to  determine  the 


16 


body  lift  and  control  surface  effects.  Model  scale  experiments  were  used  to  determine 
the  unsteady  motion  effects. 

5.1.  Experiment  Apparatus 

5.1.1.  Model  Geometry 

Figure  2  shows  the  geometry  of  the  full  scale  model.  The  small  scale  model  is 
geometrically  similar  at  a  scale  of  0.4334.  This  scale  was  selected  to  provide  the  smallest 
model  that  would  contain  the  transducer  discussed  in  Section  5.1.2  without  incidental 
contact  between  the  transducer  and  the  model. 


Ballast 


Forward 

Transducer 


Figure  2.  Full  Scale  Model  Geometry 

Figure  3  shows  the  full  scale  model  mounted  in  the  United  States  Naval  Academy 
Towing  Tank.  Figure  4  shows  the  0.4334  scale  model  mounted  in  the  Massachusetts 
Institute  of  Technology  Marine  Computation  and  Instrumentation  Laboratory. 


Figure  3.  Full  Scale  Model  Mounted  in  United  States  Naval  Academy  Towing  Tank 


Figure  4. 0.4334  Scale  Model  Mounted  at  the  MIT  Marine  Instrumentation  and  Computation 
Laboratory 

5.1.2.  Force  and  Moment  Measurement 

The  forces  and  moments  were  measured  using  a  UDW3  underwater  transducer 
manufactured  by  Advanced  Mechanical  Technology,  Inc.  The  transducer,  shown  in 
Figure  5,  is  able  to  simultaneously  measure  forces  and  moments  in  all  of  the  three 
orthogonal  directions  (making  six  measurements  of  forces  and  moments)  and  is  suitable 
for  underwater  applications.  A  pressure  compensating  bladder  in  the  transducer  equalizes 
internal  and  external  pressures  to  allow  underwater  operation  with  little  effect  of 
hydrostatic  pressure.  The  capacities  and  general  specifications  of  the  dynamometer  are 
shown  in  Table  2. 


Figure  5.  UDW3  Underwater  Sensor.8 


The  transducer  was  mounted  to  a  bulkhead  within  the  volume  of  the  vehicle.  The  strut 
was  attached  to  the  end  of  the  transducer  not  attached  to  the  vehicle.  Sufficient  clearance 
was  provided  to  ensure  the  transducer  output  was  not  compromised  by  contact  with  the 
sides  of  the  vehicle. 


Table  2.  Dynamometer  Capacity  and  Specifications9 


Vertical  and  Lateral  Force  Capacity 

556  N 

Axial  Force  Capacity 

1112.1  N 

Pitch  and  Yaw  Moment  Capacity 

28.2  N-m 

14.1  N-m 

Vertical  and  Lateral  Force  Sensitivity 

2.7  fjV  I(V  *  N) 

Axial  Force  Sensitivity 

.67  juV/(V*N) 

Pitch  and  Yaw  Moment  Sensitivity 

137.2  juV/(V*N-m) 

Roll  Moment  Sensitivity 

97.4 /jV /(V  *  N  -m) 

Vertical  and  Lateral  Force  Stiffness 

5.3  x  106  N/m 

Axial  Force  Stiffness 

7.88  x  10v  N/m 

Roll  Moment  Stiffness 

5.7  x  103  N-m/radian 

Weight 

2kg 

Recommended  Excitation 

10  V  or  less 

Crosstalk 

<  2%  on  all  channels 

Temperature  Range 

-17  to  52°  C 

Force  Channel  Hysteresis 

±  0.2%  Full  Scale  Output 

Force  Channel  Non-Linearity 

±  0.2%  Full  Scale  Output 

Excitation  for  the  transducer  and  amplification  for  the  output  were  provided  by  a 
MSA-6  Mini-Amplifier  also  developed  by  AMTI.  This  amplifier,  shown  in  Figure  6, 
provides  excitation  and  amplification  for  up  to  six  channels.  The  excitation  is  selected  by 
individual  jumpers  for  each  channel  and  ranges  from  2.5  to  10  volts.  The  gain  for  each 
channel  is  also  selectable  by  jumpers  and  ranges  from  1000  to  4000.  The  output  of  the 


19 


amplifier  is  ±10  VDC.  The  amplifier  contains  an  auto-zero  feature  that  allows  for  push 
button  zeroing  of  the  output  of  the  load  cell. 


Figure  6.  MSA-6  Mini  Amplifier10 

The  output  of  the  amplifier  was  connected  to  an  analog-to-digital  converter  installed 
in  a  notebook  computer.  The  system  control  software  sampled  the  six  channels  of  output 
of  the  load  cell  and  the  six  positions  of  the  gantry  system  at  an  operator  selected 
frequency  of  25  Hz. 

5.1.3.  United  States  Naval  Academy  Tests 

Full  scale  model  testing  was  performed  at  the  United  States  Naval  Academy 
Hydromechanics  Laboratory  shown  in  Figure  7.  This  set  of  tests  included  determining 
the  forces  and  moments  resulting  from  body  angles  in  pitch  and  yaw  and  control  surface 
angles.  The  towing  tank  used  was  120  ft  long,  8  ft  wide,  and  5  ft  deep.  The  towing  tank 
included  a  wave  making  machine,  a  wave  absorbing  beach  and  a  moving  carriage.11 


Figure  7.  United  States  Naval  Academy  Hydromechanics  Laboratory  Towing  Tank 


% 
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5.1.4.  Massachusetts  Institute  of  Technology  Tests 

Small  scale  model  testing  was  performed  at  the  Marine  Instrumentation  and 
Computation  Laboratory  at  the  Massachusetts  Institute  of  Technology.  These 
experiments  were  performed  in  order  to  determine  the  forces  and  moments  associated 
with  unsteady  motion.  The  model  was  moved  in  prescribed  sinusoidal  motions  and  the 
resultant  forces  and  moments  were  measured. 

The  laboratory  contains  a  tank  and  a  gantry  system  capable  of  simultaneous  motion  in 
five  degrees  of  freedom.  The  experimental  tank  is  10  m  long,  4  m  wide,  and  1  m  deep. 
The  gantry  system  consists  of  five  different  motors  and  several  gear  assemblies  to  ensure 
smooth  operation  at  the  speeds  and  frequencies  required  for  the  experiments.  The  gantry 
system  is  computer  controlled  for  precise  positioning. 

The  testing  and  correction  of  the  very  sophisticated  gantry  system  and  control 
software  occupied  a  significant  amount  of  the  time  allocated  for  the  performance  of  this 
research.  The  gantry  and  control  software  was  designed  and  assembled  by  D’Ambra 
Technologies,  the  only  firm  known  to  the  research  supervisor  to  be  capable  of  developing 
the  system  and  the  software.  The  original  contract  called  for  completion  of  the  gantry 
system  by  January  2002. 

Testing  began  in  June  2002  and  significant  problems  with  the  system  and  Control 
software  were  soon  identified.  Correction  of  the  problems  introduced  several  weeks  of 
delay.  The  cycle  of  problem  identification  and  correction  continued  until  very  early  in 
2003.  In  this  process  the  vertical  axis  controls  were  completely  redesigned.  The  original 
stepper  motors  were  found  to  be  inadequate  and  were  replaced  by  servo  motors.  The 
pitch  mechanism  was  strengthened  three  times  to  be  able  to  provide  the  desired  frequency 
and  amplitudes  of  oscillation.  The  gantry  system  and  control  software  were  believed  to  be 
reliable  and  accurate  in  early  April  2003. 

The  research  team  also  encountered  problems  related  to  unidentified  faults  in  the  force 
measurement  system.  AMTI  conducted  significant  troubleshooting  of  the  load  cell  and 
connections  on  several  occasions  to  determine  the  cause  of  the  abnormal  readings.  Some 
of  the  abnormal  readings  were  attributed  to  a  fault  in  the  cabling  and  others  were 
attributed  to  the  high  level  of  electronic  noise  in  the  long  cables  of  the  system. 

Eventually,  all  of  the  issues  with  the  force  measurement  system  were  corrected. 

One  difficult  issue  with  the  force  measurement  system  that  was  discovered  late  in  the 
process  is  that  the  measured  forces  and  moments  often  represent  a  very  small  fraction  of 
the  capability  of  the  transducer.  The  manufacturer  states  that  the  transducer  provides 
accurate  results  at  very  small  fractions  of  its  capability,  but  this  needs  experimental 
verification. 

Several  hundred  experiments  were  conducted  during  the  process  of  identifying  and 
correcting  system  problems.  Experiments  were  performed  by  the  author,  by  an 
independent  contractor,  and  by  several  undergraduate  students  acting  with  supervision. 
The  data  from  these  experiments  needs  to  be  closely  examined  to  determine  if  the 
experiments  are  valid.  The  first  future  work  that  will  be  done  is  to  predict  what  those 
tests  should  show  and  check  for  agreement.  If  the  data  are  found  to  be  valid,  they  will 
contribute  to  a  more  complete  data  set  with  less  variance  that  will  allow  for  better 
modeling  of  the  vehicle  behavior. 
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5.2.  Design  of  Experiments 

The  variables  that  affect  the  behavior  of  an  underwater  vessel  include  the  depth  of  the 
water,  the  submergence  of  the  vessel,  the  forward  velocity,  and  the  frequency  and 
amplitude  of  oscillation.  The  Central  Composite  Method  was  used  for  Design  of 
Experiments  in  order  to  reduce  the  total  number  of  experiments  required. 

The  Central  Composite,  or  Box-Wilson,  Design  is  a  three-  or  five-level  design  that 
includes  the  comer,  center,  and  axial  points  of  the  design  space.  The  three-factor  Central 
Composite  Design  space  is  shown  in  Figure  8. 


Velocity 


Figure  8.  Central  Composite  Method 


The  three  factor  design  space  is  developed  from  1 5  point  designs:  a  center  point 
design,  eight  comer  point  designs,  and  6  axial  point  designs.  This  model  represents  the 
response  surface  more  accurately  than  most  other  methods  since  the  comer  points  are 
included.  Comer  points  represent  the  limits  of  the  experimental  space.  However, 
attempting  to  reach  these  comer  point  designs  may  strain  the  engineering  model12. 

The  selection  of  test  points  in  an  incomplete  matrix  on  the  basis  of  orthogonal  numeric 
functions  is  fine  when  the  dependent  variable  depends  linearly  on  the  input  variables. 
However,  for  things  like  a  nonlinear  relation  between  force  coefficient  and  excitation 
frequency,  it  is  better  to  be  sure  that  all  comers  in  the  test  space  are  tested  so  that  the 
mathematical  model  will  interpolate  rather  than  extrapolate. 

Several  experiments  that  had  been  planned  were  not  performed  due  to  limitation  of  the 
gantry  system  at  higher  speeds  and  higher  frequencies  of  oscillation.  Other  experiments 
were  not  performed  due  to  physical  constraints  of  the  gantry  system.  Appendix  A  lists 
the  full  scale  experiments  that  were  performed.  Appendix  B  lists  the  model  scale 
experiments  that  were  used  for  analysis.  All  submergences  listed  in  the  test  matrices  are 
to  the  center  of  the  body. 

A  large  number  of  other  experiments  were  performed  earlier  in  the  test  program,  but 
uncertainty  in  the  equipment  behavior  resulted  in  uncertainty  in  the  quality  of  the  data 
and  data  from  those  experiments  was  not  used.  That  data  will  be  reevaluated  as  part  of 
future  work. 

5.3.  Analysis  of  Experimental  Data  to  Extract  Measured  Forces  and 
Moments 

For  testing  performed  at  MIT,  the  conversion  from  raw  force,  moment,  and  position 
data  was  performed  using  MATLAB  routines  developed  by  the  author.  The  transducer 
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provided  data  on  forces  and  moments  in  the  form  of  voltages  for  each  channel  that  had  to 
be  converted  to  the  MKS  system.  The  gantry  system  provided  data  on  the  position  of  the 
system  in  a  numerical  format  that  had  to  be  converted  to  the  MKS  system  for  analysis. 

For  testing  performed  at  the  United  States  Naval  Academy,  raw  force  and  moment 
data  were  converted  using  routines  developed  the  author  and  by  LTJG  Greg  Sabra, 

USCG.  The  steady  force  results  were  determined  in  a  manner  very  similar  to  the  method 
used  to  compute  steady  force  test  results  for  MIT  tests.  The  method  for  computing  MIT 
test  results  will  be  discussed  in  a  later  section.  A  complete  discussion  of  LTJG  Sabra’s 
code  is  contained  in  his  thesis  entitled  ““Wave  Effects  on  Underwater  Vehicles  in 
Shallow  Water.” 

All  of  the  MIT  test  conditions  were  listed  in  a  common  Excel  file  called  “MIT  Test 
Plan.xls”.  This  file  contains  separate  worksheets  for  each  of  the  many  series  of  tests  that 
were  performed.  These  worksheets  look  very  similar  to  the  table  contained  in  Appendix 
B,  with  the  addition  of  two  columns  at  the  left  to  record  date  and  time  information  for 
each  experiment.  The  worksheets  in  the  MIT  Test  Plan  file  were  used  by  the  MATLAB 
routines  to  determine  which  data  files  to  analyze  and  what  some  of  the  test  conditions 
were  for  each  experiment.  The  test  conditions  obtained  from  the  test  plan  were  the  water 
depth  and  submergence  of  the  vehicle  during  the  test.  All  other  test  conditions  were 
extracted  directly  from  the  test  data  file. 

5.3.1.  User  Interface  and  Data  File  Management 

The  user  interface  and  the  file  management  were  performed  by  a  MATLAB  routine 
called  “AutoanalyzeXls.m”.  This  file  is  contained  as  Appendix  C.  After  the  user  starts 
this  program,  the  user  selects  the  series  of  experiments  to  be  analyzed  by  entering  the 
number  corresponding  to  the  desired  series.  All  series  that  have  been  performed  are 
listed,  even  those  that  are  suspected  to  be  of  little  value.  After  the  test  series  has  been 
selected,  the  program  imports  the  list  of  experiments  and  the  depth  and  submergence 
information.  The  program  then  calls  other  MATLAB  routines  to  analyze  the  data  files. 
For  steady  force  tests,  the  analysis  program  is  "AnalyzemodXlsSF.m".  For  all  other  tests, 
the  analysis  program  is  "AnalyzemodXls.m". 

The  analyses  were  performed  using  the  data  files  recorded  by  the  notebook  computer 
in  Excel  format.  An  example  of  a  data  file  is  included  as  Appendix  D. 

5.3.2.  Steady  Force  Data  Analysis 

Steady  force  tests  at  MIT  involved  towing  the  vehicle  down  the  tank  with  a  steady 
angle  of  yaw  or  pitch.  These  tests  were  analyzed  using  "AnalyzemodXlsSF.m", 
contained  in  Appendix  E.  "AnalyzemodXlsSF.m"  starts  by  importing  the  data  file 
identified  by  “AutoanalyzeXls.m”.  The  program  determines  the  ordered  parameters  and 
the  date  and  time  at  which  the  experiment  occurred.  Then,  the  file  eliminates  the  first  1 .2 
seconds  of  data  to  allow  for  gantry  acceleration  and  any  data  recorded  after  the  gantry 
velocity  returns  to  zero  at  the  end  of  the  test.  The  remaining  position  data  is  converted  to 
the  MKS  system  using  a  conversion  factor.  The  remaining  force  and  moment  data 
undergo  a  more  detailed  analysis. 

The  voltage  output  of  the  transducer  is  converted  to  forces  and  moments  using 

F  = - ^ - r  (14) 

Gain  x  Vexc  x  S'  x  1CT6 
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where 


F  is  the  calculated  force  or  moment 
Vout  is  the  output  voltage  recorded  by  the  computer 
Gain  is  the  gain  of  that  channel  in  the  amplifier 
Vcxc  is  the  excitation  voltage  of  the  channel,  and 
S  is  the  sensitivity  of  the  channel. 

The  calculation  of  equation  (14)  is  performed  using  matrices  so  that  the  effects  of  cross¬ 
talk  in  the  transducer  can  be  accounted  for. 

The  mean  force  is  calculated  by  taking  the  mean  of  the  forces  measured  in  the  data 
interval  and  shifting  the  origin  of  the  mean  force  from  the  origin  of  the  load  cell  to  the 
origin  of  the  vessel  coordinate  system,  defined  to  be  at  the  midships  of  the  vessel.  The 
origin  shift  was  done  using 


^ midships  ^  transducer  transducer 

N  =  N  +Yx 

midships  transducer  transducer 


(15) 


The  origin-shifted  mean  forces  and  moments  were  written  to  a  common  output  file  that 
contained  the  mean  force  and  moment  data  for  all  analyzed  experiments.  This  output  file 
and  explanatory  notes  are  included  as  Appendix  F. 


5.3.3.  Analysis  of  Experiments  Involving  Unsteady  Motion 

The  analysis  of  experiments  involving  unsteady  motion  was  performed  using 
“AnalyzemodXls.m”  called  by  “AutoanalyzeXls.m”.  The  code  is  included  as  Appendix 
G.  This  is  the  most  complicated  of  the  codes  used  for  this  research  and,  as  a  result,  is  the 
most  heavily  commented. 

The  first  section  of  the  code  identifies  and  defines  most  of  the  variables  used  in  the 
code.  Next,  the  code  initializes  by  reading  the  data  file  and  gathering  some  basic 
information  about  the  parameters  of  the  experiment.  The  last  four  lines  of  the  data  file 
contain  information  about  the  ordered  frequency  and  amplitudes  of  oscillation  as  well  as 
sample  frequency,  velocity,  and  travel  duration  and  distance.  Then,  the  actual  sample 
frequency  is  calculated  by  taking  the  inverse  of  the  average  interval  between  data  points 
according  to 

f sample  /  . 

mean(At) 

The  ordered  sample  frequency  was  always  25  Hz,  but  for  a  certain  period  of  time  during 
the  research  errors  in  the  control  software  resulted  in  data  being  taken  at  other 
frequencies. 

The  code  drops  the  first  1 .2  seconds  of  data  to  allow  for  the  acceleration  of  the  gantry. 
The  time  interval  to  drop  was  chosen  short  enough  to  allow  sufficient  time  remaining  to 
have  several  periods  of  oscillation  remaining  but  long  enough  to  remove  the  majority  of 
the  acceleration  period.  The  code  also  drops  data  recorded  after  the  vessel  completed  its 
travel  along  the  tank.  This  was  necessary  because  the  control  software  continued  to 
collect  data  until  the  ordered  time  period  of  the  experiment  was  completed,  whether  or 
not  the  travel  distance  had  been  accomplished.  Setting  the  time  period  of  travel  too  short 
resulted  in  sudden  stops  of  the  gantry  causing  large  accelerations  on  both  the  vehicle  and 
the  gantry  system.  Ordered  durations  were  made  longer  than  absolutely  necessary  to 
prevent  this  mechanical  shock  to  the  system  and  prolong  the  life  of  the  apparatus. 


The  analysis  code  determines  the  frequency  of  oscillation  by  checking  the  input 
parameters  in  the  data  file  to  determine  the  ordered  frequency  of  oscillation.  This 
information  is  used  to  ensure  that  the  data  to  be  analyzed  consisted  of  an  integer  number 
of  wavelengths  of  the  oscillation.  This  feature  was  absolutely  necessary  to  get  highly 
accurate  results  from  the  Fourier  analysis  that  takes  place  later  in  the  program.  The 
period  of  a  cycle  is  given  by 

period  = - - - 

frequency 

The  duration  of  data  recorded  was  found  by  taking  the  difference  in  time  between  the 
first  and  last  remaining  data  points.  The  number  of  periods  recorded  is 

,  „  .  ,  duration 

number  of  penods  = - - 

period 

The  number  of  data  points  retained  for  Fourier  analysis  is  found  by  rounding  down  to  the 
next  integer  the  product  of  sample  frequency,  period,  and  the  number  of  periods 
according  to 

#  of  data  points  =  round (fsampIe  *  period  *  floor{ number  of  periods)) 

“round”  is  a  MATLAB  function  that  rounds  the  element  to  the  nearest  integer,  “floor”  is 
a  MATLAB  function  that  round  the  element  to  the  next  lower  integer. 

Once  an  integer  number  of  data  points  is  established  it  the  mean  force  and  moments 
are  subtracted  from  all  measured  forces  and  moments  in  order  to  remove  steady  effects. 

Next,  the  voltages  from  the  transducer  are  converted  to  forces  and  moments  using 
equation  (14)  and  the  numeric  position  data  is  converted  to  metric  system  position  data 
using  known  relationships  between  the  controller  data  and  gantry  motion. 

The  force,  moment  and  location  data  are  conditioned  by  the  program  in  preparation 
for  the  Fourier  analysis.  With  the  position  data  in  metric  format,  the  program  translates 
the  position  data  from  the  location  of  the  strut  to  the  vehicle  midships.  For  linear  motion, 
the  motion  of  the  strut  forward  of  midships  represents  the  motion  of  midships.  For 
angular  motion,  this  is  not  the  case.  The  effect  of  angular  motion  on  the  x,y,and  z 
position  of  midships  is  calculated  by 

y midships  -  y strut  -  distance  sin 


where 


^ midships  ^ strut  ^ 'strut  ^  ^ 


e 


n 

180 


•7  —  Z  _  J 

^ midships  ~  strut  ^ strut 


1  -  cos 


V 


(17) 


6  is  the  pitch  angle 
y/  is  the  yaw  angle 

distance  is  the  distance  along  the  x  axis  from  the  strut  to  midships. 

Lstrut  is  the  length  of  the  strut  arm  from  its  pivot  point  to  the  vehicle. 

The  forces  and  moments  are  shifted  from  having  their  origin  at  the  transducer  to  having 
their  origin  at  midships  using  equation  (15).  Also,  the  data  is  interpolated  into  even 
intervals  of  exactly  0.4  seconds.  The  mean  value  of  position  for  each  channel  except  the 


x  position  is  subtracted  to  remove  any  bias  in  the  position  data.  Then,  the  data  matrix  is 
padded  with  zeros  to  obtain  exactly  2048  data  points. 

The  most  precise  Fourier  transformation  requires  the  data  to  have  an  integer  number 
of  periods  of  the  waveform  and  the  frequency  of  signal  to  be  analyzed  must  be  a  multiple 
of  the  fundamental  frequency  of  the  data  sample.  For  an  interval  of  2048  data  points 
being  sampled  at  25  Hz,  the  fundamental  frequency  is 

nu^poi^  =  i^  =  ft012207  <18) 

The  frequencies  of  oscillation  used  for  this  research  are  0.402831, 0.79346,  and  1.19629 
Hz  representing  33,  65  and  98  times  the  fundamental  frequency  of  the  analysis.  The 
sample  frequency  is  assured  by  interpolating  the  data  into  exact  time  intervals  between 
data  points.  As  a  result,  the  force,  moment,  and  position  data  relating  to  oscillation  are 
readily  extracted  using  Fourier  analysis. 

The  first  step  in  the  Fourier  analysis  is  to  begin  to  build  the  data  matrix  by 
constructing  the  frequency  column.  The  first  row  is  assigned  a  frequency  of  zero  Hz  and 
each  successive  row  is  assigned  a  frequency  of  the  row  number  multiplied  by  the 
fundamental  frequency.  The  fast  Fourier  transformation  is  applied  to  the  force,  moment, 
and  position  data.  To  compensate  for  the  zero  padding  added  earlier,  the  value  of  each  of 
the  Fourier  coefficients  is  multiplied  by  the  ratio  of  the  padded  size  to  the  unpadded  size. 

After  the  Fourier  transformation  occurs,  each  coefficient  has  both  a  magnitude  and 
phase  associated  with  it.  The  phase  of  each  of  the  coefficients  is  changed  to  make  it 
relative  to  the  phase  of  the  motion  that  produced  the  force.  This  is  done  by  multiplying 
every  coefficient  by  e~l<p  where  <p  is  the  phase  angle  of  the  driving  motion  at  that 
frequency. 

Low  pass  filters  were  installed  on  all  of  the  force  data  collection  channels  to  reduce 
the  effects  of  electronic  noise  in  the  system.  The  effects  of  these  filters  is  removed  from 
each  channel  of  force  and  moment  using 

1  unfiltered  =  V filtered6  ^  0^) 

where  T}unfiUered  is  the  amplitude  of  the  signal  with  filtering  effects  removed 

rt filtered  *s  the  amplitude  of  the  signal  after  filtering 

(o  is  the  frequency  of  oscillation 
C  is  the  capacitance  of  the  filter,  and 
R  is  the  resistance  of  the  filter. 

The  effects  of  filtering  at  the  frequencies  of  oscillation  altered  the  magnitude  by  less  than 
one  percent  and  the  phase  angle  by  less  than  five  degrees. 

The  program  determines  the  frequency  of  motion  closest  to  the  ordered  frequency  of 
oscillation  for  all  motions.  The  actual  amplitudes  of  motion  and  forces  are  converted 
back  into  the  time  domain  from  the  frequency  domain  using 


7(0  = . 


7(<y) 


where  77(f)  is  the  amplitude  in  the  time  domain 

rj(co)  is  the  amplitude  in  the  frequency  domain,  and 
2048  is  the  number  of  data  points  used  in  the  analysis. 
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The  actual  amplitudes  of  motion  are  used  to  calculate  the  inertial  forces  experienced 
by  the  vehicle  using  equation  (2).  The  inertial  forces  are  then  subtracted  from  the 
measured  forces  and  moments  to  leave  only  the  hydrostatic  forces  and  moments 
remaining.  The  mass,  inertia,  and  center  of  gravity  terms  used  to  calculate  the  inertial 
forces  and  moments  were  derived  by  performing  oscillation  tests  in  air.  The  model  was 
filled  with  water  to  ensure  that  the  effective  mass  of  the  vehicle  was  the  same  as  it  would 
be  if  the  vehicle  were  in  water.  This  guaranteed  that  the  only  significant  forces  measured 
were  the  inertial  forces.  The  validity  of  the  process  was  checked  by  performing 
oscillation  tests  in  air  and  subtracting  the  calculated  inertial  components.  The  results  of 
these  inertial  calculation  checks  are  included  as  Appendix  H. 

At  this  point  in  the  program  all  of  the  hydrodynamic  forces  and  moments  at  the 
frequencies  of  oscillation  and  their  phases  relative  to  the  driving  motion  are  determined. 
The  results  are  conditioned  in  order  to  have  positive  amplitudes  and  have  magnitudes  of 
the  phase  angles  less  than  180  degrees. 

The  program  has  also  determined  the  forces  and  moments  at  twice  and  three  times  the 
oscillation  frequency.  In  all  cases  the  forces  and  moments  at  multiples  of  the  oscillation 
frequency  have  magnitudes  of  approximately  10%  or  less  of  the  forces  and  moments  at 
the  oscillation  frequency.  This  indicated  the  response  of  the  complete  hydrodynamic 
system  is  linear  and  that  non-linear  forces  and  moments  can  safely  be  neglected  in 
analyzing  vehicle  dynamics  or  in  designing  control  systems.  This  also  indicates  that 
problems  such  as  hysteresis  in  the  load  cell  o-rings  were  unlikely. 

The  results  are  written  to  an  output  file  and  contain  all  of  the  test  information  as  well 
as  the  actual  amplitudes  and  frequencies  of  oscillation  and  the  amplitudes  and  frequencies 
of  all  six  forces  and  moments.  The  phase  angle  between  the  force  or  motion  and  the 
driving  is  also  listed.  The  output  also  includes  the  frequency,  amplitude,  and  phase 
information  for  the  second  and  third  harmonics  forces  and  moments.  A  partial  example 
of  the  output  file  is  contained  in  Appendix  I.  A  complete  output  file  consists  of  one  row 
of  103  columns  for  each  file  analyzed. 

5.3.4.  Curve  Fitting  the  Experimental  Results 

The  following  sections  describe  how  the  force  and  moment  results  were  analyzed  and 
present  the  resulting  maneuvering  coefficients.  The  effects  of  submergence  and  speed  on 
the  maneuvering  coefficients  will  be  analyzed.  Part  of  that  analysis  will  involve  curve 
fitting  the  data  to  determine  the  functional  relationships  involved.  The  program  used  to 
do  the  curve  fitting  is  a  MATLAB  routine  generated  by  the  author  called 
“CoeffSolver.m”.  The  code  is  included  as  Appendix  J. 

“CoeffSolver.m”  uses  user-coded  values  of  the  coefficients  and  the  test  parameters  of 
Length/Submergence  and  Froude  Number  to  perform  a  least  squares  regression  of  the 
data.  When  sufficient  test  data  is  present  a  second  order  equation  is  derived.  When  there 
is  not  sufficient  data  for  a  second  order  equation  a  first  order  equation  is  used. 

In  general,  the  set  of  equations  is  of  the  form 

a,,*,  +  ai2x2  +  anx2  =  bx 

a2lXl  ^22X2  ®23X3  =  ^2 

a31Xl  ®32X2  a33X3  =  ^ 3 

which  can  be  written  in  matrix  form 


Ax-B 

where 

A  is  an  m  x  n  matrix 
x  is  an  array  of  size  m 
B  is  an  array  of  size  m. 

Matrix  A  represents  the  known  parameters  of  the  equations  such  as  the  speed  and 

submergence  of  the  vehicle  for  the  experiment.  Also,  the  terms  an  usually  equal  1  to 

allow  for  some  constant  to  be  built  into  the  equation.  The  B  array  represents  the 

measured  or  calculated  data  points.  The  array  x  represents  the  coefficients  that  are 

calculated  to  best  represent  the  data  using  a  least  squares  linear  regression.  The  least 

squares  regression  calculates  the  elements  of  array  x  that  will  minimize  the  sum  of  the 

squares  of  the  errors  between  the  predicted  and  the  calculated  values.  This  manipulation 

is  performed  easily  in  MATLAB  using  x  =  A\B  . 

To  measure  the  quality  of  fit,  the  program  finds  the  root  mean  square  of  the  difference 

between  the  predicted  and  the  measured  coefficients.  The  code  uses 

_  „  (Pr  edicted  -  Actual) 

Difference  = - 

Actual 

and 


Difference rms 


norm  {Differ en ce ) 
sqrt{ #  of  elements) 


where  the  norm  of  the  Difference  array  is  the  largest  singular  value  in  the  Difference 


array. 

The  output  of  the  Coefficient  Solver  Program  is  included  as 


5.4.  Analysis  of  Experimental  Results 


5.4.1.  Rudder  and  Stern  Planes  Effects 

The  effects  of  the  rudder  and  stem  planes  were  examined  by  performing  tests  at  the 
United  States  Naval  Academy  with  the  control  surfaces  at  no  angle  and  with  the  control 
surfaces  deflected  to  seven  degrees.  The  results  were  normalized  by  subtracting  the  lift 
and  moment  at  zero  degrees  from  the  lift  and  moment  measured  with  control  surface 
deflection  to  remove  any  imbalance  in  loading.  The  lift  coefficient  per  degree  of  fin 
deflection  was  calculated  according  to 

c, v-tA*—  PD 

PU2ArJ 

where  Cl8  is  the  lift  coefficient  per  degree,  Lift  is  the  measured  lift  force,  U  is  the  vehicle 
forward  velocity,  Afm  is  the  effective  area  of  the  rudder,  and  8  is  the  control  surface 
angle.  In  order  to  provide  more  data  points  from  the  same  experiments,  the  measured 
moments  were  converted  into  pseudo-lift  forces  by  dividing  the  moments  by  the  length  of 
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the  control  surface  moment  arm,  0.7  m.  These  were  then  converted  into  coefficients 
using  equation  (21). 

For  high  aspect  ratio  wings,  with  aspect  ratios  greater  than  five,  the  theoretical 
approximation  found  by  Hoemer13  is 


=  dC L_ 

da 


10  + 


20 
yf  Re 


(22) 


where  ARC  is  the  aspect  ratio  found  by 


A  Re  = 


r  Span2"' 
Area 


(23) 


The  area  was  estimated  by  combining  the  calculated  area  of  each  of  the  fins  with  the 
estimated  effective  area  provided  by  the  body  between  the  fins.  Figure  9  shows  the 
control  surface  geometry.  Using  an  area  of  0.02  m2  and  span  of  0.254  m,  the  aspect  ratio 
is  3.23,  and  Cl  is  0.0618/Degree. 


Figure  9.  Control  Surface  Geometry 


For  aspect  ratios  between  three  and  five  Hoemer  recommends  using 


■'La 


da 


10  + 


10 
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(24) 


A  R2  yfR 

which  yields  Cl  of  0.053/Degree.14 

Figure  10  shows  the  measured  rudder  lift  coefficients  plotted  against  the  ratio  of  body 
length  to  submergence.  The  figure  also  shows  the  linear  approximation  to  the  data  and 
the  theoretical  value  for  intermediate  aspect  ratio  fins  in  an  infinite  fluid.  The  data  for 
Figure  10  are  presented  in  Table  3. 

The  linear  approximation  to  the  data  is 


=  -0.00 1 3  +  0.0518 


Submergence 


(25) 


The  value  of  the  linear  approximation  at  deep  submergence  is  0.05 1 8/Degree,  very  near 
the  theoretical  value  of  0.053/Degree  calculated  using  equation  (24). 
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Figure  10.  Rudder  Lift  Coefficient  for  Various  Values  of  Length/Submergence 


Table  3.  Rudder  Lift  Coefficient  Data 


Test 

Speed 

L/Subm. 

Rudder 

Angle 

Rudder 

Lift 

Normalized 
Rudder  Lift 

Rudder 
Lift  Coeff 

Yaw 

Moment 

Normalized 

Yaw 

Moment 

Yaw 

Moment 

Coeff 

m/s 

Degrees 

N 

N 

1 /Degree 

N-m 

N-m 

1 /Degree 

Measured  Lift 

N1 

2.06 

2.07 

0 

-0.915 

0.000 

0.333 

0 

N6 

0.515 

4.14 

0 

-0.259 

0.000 

0.302 

0.000 

N7 

1.03 

4.14 

0 

-1.729 

0.000 

1.336 

0.000 

N8 

2.06 

4.14 

0 

-7.349 

0.000 

6.402 

0.000 

N9 

2.06 

4.14 

7 

12.267 

19.616 

0.0662 

-17.316 

0.0825 

N15 

0.515 

8.26 

0 

0.280 

0.000 

0.196 

0.000 

N16 

2.06 

8.26 

0 

2.810 

0.000 

1.512 

0.000 

N22 

2.06 

8.26 

7 

8.001 

5.191 

0.0175 

-10.594 

-12.106 

0.0577 

2.06 

8.26 

7 

7.922 

5.112 

0.0172 

BEB 

-14.595 

0.0695 

N9m 

2.06 

4.14 

7 

-7.832 

0.0264 

Lift  Calculated 
from  Moment 

N22m 

2.06 

8.26 

7 

-17.294 

0.0583 

N24m 

2.06 

8.26 

7 

-20.850 

0.0703 

The  effects  of  stem  planes  were  determined  in  a  manner  similar  to  those  of  the  rudder 
with  similar  results.  Figure  1 1  shows  the  measured  Stem  Planes  Lift  Coefficients  plotted 
against  the  ratio  of  body  length  to  submergence.  Figure  1 1  also  shows  the  linear 
approximation  to  the  data  and  the  calculated  theoretical  value.  The  data  for  Figure  1 1  are 
included  as  Table  4. 
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Figure  11.  Stern  Planes  Lift  Coefficient  for  Various  Values  of  Length/Submergence 


Table  4.  Stern  Planes  Lift  Coefficient  Data 


Test 

Speed 

L/Subm. 

Stem 

Planes 

Angle 

Stern 

Planes 

Lift 

Normalized 
Stem 
Planes  Lift 

Stern 
Planes 
Lift  Coeff 

Pitch 

Moment 

Normalized 

Pitch 

Moment 

Pitch 

Moment 

Coeff 

m/s 

N 

1 /Degree 

N-m 

N-m 

1/Degree 

Measured  Lift 

an 

0 

2.850 

0.000 

-3.227 

H 

whesm 

7 

19.131 

16.281 

0.0549 

4.982 

0.0391 

mt 

PI  B 

4.14 

0 

n 

0.000 

-0.157 

H 

4.14 

0 

mm 

-0.778 

fcnsfl 

0 

-6.761 

N10 

| 

-7 

n 

0.0601 

-8.460 

-1.699 

0.0081 

sa 

0 

sREm 

-0.162 

an 

Btl!l 

0 

fin 

-8.868 

Em 

-7 

0.0510 

RIBS! 

PB 

n 

-7 

| 

0.0712 

HH 

li  Vi  Will 

Lift  Calculated 
from  Moment 

m 

2.06 

2.07 

7 

0.0396 

N20m 

2.06 

8.26 

-7 

0.0563 

N24m 

2.06 

8.26 

-7 

0.0712 

These  results  seem  to  show  that  the  rudder  has  a  slightly  reduced  effect  when  near  the 
surface  and  the  stem  planes  have  a  slightly  greater  effect  when  near  the  surface.  In  both 
cases,  the  linear  approximation  of  the  data  closely  approaches  the  theoretical  value  for  a 
submerged  body  in  an  infinite  fluid.  The  opposite  slopes  of  the  linear  approximations  are 
not  explainable  at  this  time.  The  small  variations  shown  in  the  coefficients  over  the  range 
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of  submergences  monitored  show  that  the  effect  of  submergence  is  very  small  for 
submergences  greater  than  10%  of  body  length. 


5.4.2.  Sway  Force  and  Yaw  Moment  Due  to  Sway  Motion 

Experiments  to  determine  the  effects  of  sway  motion  were  performed  in  various 
combinations  of  submergence,  velocity,  frequency  and  amplitude  of  oscillation.  Table  5 
contains  the  test  conditions  and  key  results. 


Table  5.  Test  Conditions  and  Results  for  Sway  Force  and  Yaw  Moment  Due  to  Sway  Motion 


Sway  Motion 

Hydrodynamic  Sway  Force,  Y  j 

Frequency 

Amplitude 

■aimlE 

m&XM 

Yvdot 

Yvdot’ 

USSSSSM 

Y„ 

Yv’ 

Hz 

m 

N 

N 

_ kS _ 

N 

WE mm 

0.543 

0.40283 

0.1 

2.217 

-0.0175 

-1.1944 

-4.7044 

0.252 

Mil 

0.40283 

0.1 

2.641 

mi*  wtM 

-1.6937 

0.398 

MiTMM 

0.79346 

0.1 

8.644 

WiSHl 

-2.7725 

-5.0439 

gillie  ■ 

0.543 

■AsHdi 

0.40283 

0.1 

2.677 

-38.1 

-3.2812 

-1.6517 

giliW 

0.252 

0.40283 

0.1 

2.216 

im 

-2.9156 

-0.0175 

-1.1841 

-4.6652 

WEEiMM 

I  Sway  Motion 

Hydrodynamic  Yaw  Moment,  N  1 

Frequency 

Phase 

WE3SM 

Nvdot 

Nvdot’ 

■C9H 

Nv’ 

m 

m/s 

Hz 

m 

N-m 

■D8SH 

_ _ 

■SEH 

■Si 

0.543 

0.40283 

0.1 

0.085 

-45.4 

EHI 

EEMiliil 
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■Q HI 

0.40283 

0.1 

0.351 

0.0513 

E&U&H 

-0.34714 

-1.36874 

B.M.iWtl 

0.398 

0.79346 

0.1 

0.352 

-55.8 

msnsnm 

BiliTiliKEl 

Bililignai 

0.543 

0.40283 

0.1 

0.371 

-80.1 

■iXtlsKUM 

■****■ 

-1.44166 

0.252 

0.40283 

0.1 

0.085 

1» 

E'l'MJU 

EE®a 

1«EI 

The  hydrodynamic  forces  and  moments  are  determined  by  subtracting  the  inertial 
forces  and  moments  from  the  measured  forces  and  moments  for  each  experiment.  For 
pure  sway  motion,  the  applicable  equations  of  motion  from  equation  (9)  are 

Y  =  Yy  +  Y,v 


(26) 


(27) 


N  =  Nvv  +  Niv 

The  notation  of  complex  equations  is  used  to  separate  the  components  of  the  measured 
force  and  moment  into  the  component  related  to  velocity,  a  damping  component,  and  the 
component  related  to  acceleration,  an  added  mass  component.  This  is  done  using 

Y  =  Re(7)  +  i  Im(F)  =  Y  cos(^)  +  iY  sin(^) 

N  =  Re(  N)  +  i  Im  (TV)  =  N  cos(^)  +  iN  sin(^) 

and  understanding  that 

Re(7)  =  Y  cos(^)  =  Y^v 
Im(7)  =  Y  sin(^)  =  yvv 
Re(N)  =  Ncos(0)  =  Y-v 
Im(7V)  =  N  sin(^)  =  Yvv 


(28) 


^  where  <j>  is  the  phase  angle  taken  by  subtracting  the  phase  angle  of  the  force  from  the 

phase  angle  of  the  driving  motion. 

Sinusoidal  motions  can  be  described  by 

77  =  rjei0)t  (29) 

where  7  is  a  time  varying  position 

#  77  is  the  amplitude  of  the  sinusoidal  oscillation,  and 

c 0  is  the  frequency  of  oscillation. 
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The  velocity  component  is 


(30) 


—  =  coAe™ 
dt 


and  the  acceleration  component  is 


d2rj 

It2 


=  -CD  Aeu 


(31) 


In  the  standard  method  of  notation  the  exponential  term  is  understood  and  velocity  is 
represented  by  coA  and  the  acceleration  term  is  represented  by  -co2  A . 

Therefore, 

v  =  co  A 


(32) 


V  =  -co2  A 

By  combining  equations  (28)  and  (32),  we  obtain  the  necessary  equations  to  compute  the 
hydrodynamic  forces  and  moments  due  to  sway  according  to  . 

Im(7) 


n 

n  = 

N,= 


coA 

Re(7) 

-CD2  A 
Im(iV) 
coA 

R  e(AQ 
-co1  A 


(33) 


The  results  are  non-dimensionalized  according  to  the  method  of  Section  4. 

5.4.3.  Heave  Force  and  Pitch  Moment  Due  to  Heave  Motion 

Experiments  to  determine  the  effects  of  heave  motion  were  also  performed  in  various 
combinations  of  submergence,  velocity,  frequency  and  amplitude  of  oscillation.  Table  6 
contains  the  test  conditions  and  key  results. 

Table  6.  Test  Conditions  and  Results  for  Heave  Force  and  Pitch  Moment  Due  to  Heave  Motion 


|  Heave  Motion 

Hydrodynamic  Heave  Force,  Z  | 

BSSSB9 

EStEBl 

Amplitude 
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Zwdot' 

lm(Z) 

Zw 

Zw' 

m 
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Hz 

m 

N 

N 

_ kg 

N 

WEETM 
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■flkkkM 

0.40283 

0.1 
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WExmm 

giliTSEEl 
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Hi:  till 

0.40283 
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j 

WKSSEEMM 

-3.8243 

-0.0230 

-6.9812 

giliWtll 

0.398 

0.79346 

0.1 

6.9747 

-2.7856 

-0.0167 

EEESfl 

E>K#1 

0.488 

■EMI 

0.40283 

0.1 

2.687 

WKEE3B&M 

-3.4543 

-1.5100 

-5.9418 

EE3H 1 

0.272 

0.40283 

0.1 

2.633 

>28.9 

2.3054 

■SEE* 

-1.2727 

Efgfill 

gililiCTl 

Heave  Motion 

Hydrodynamic  Pitch  Moment,  M  1 

Submergence 

BM5SBIHI 

ESMSa 

■E3C3M 

Mwdot 

Mwdot’ 

ekeb 

Mw’ 

m 

m/s 

Hz 

m 

iiSBffSHI 

N-m 

*9 

■2ESH 

KB 

heejssh 

0.40283 

Boom 

■MkEkiM 

mmsm 

0.32703 

mmmm 

1.02217 

0.272 

0.40283 

0.1 

0.395 

■mETilil 

-0.1472 

0.23097 

■ililiVJtlil 

1.45422 

IiIiIiliMI 

0.398 

0.79346 

0.1 

0.966 

■ESI 

-0.8223 

0.32843 

1.00733 

EEB3 

0.488 

■I«W 

0.40283 

0.1 

gfigfcM 

mkesm 

■E3EZ323I 

EB 

I&bIHil 

ees sa 

1.75554 

0.272 

0.333 

0.40283 

0.1 

0.293 

BHI 

-0.2372 

0.37125 

0.00322 

0.17110 

imiKwci 

The  analysis  of  the  results  in  heave  motion  follows  t 


he  same  train  of  reasoning  as 


described  previously  for  sway  motion.  The  equations  for  heave  motion  become 


Z  =  Zvvw  +  Z> 

M  =  Mww  +  M^w 


(34) 
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and 


Im(Z) 


7  Rg(g) 
*  -co2A 

M 

w 

M 

^  2  a 

-co  A 


(35) 


5.4.4.  Discussion  of  Heave  and  Sway  Motion  Results 

The  results  of  these  experiments  are  in  agreement  with  the  expected  results.  Because 
the  sway  and  heave  motions  with  forward  velocity  create  effective  angles  of  attack  of  the 
body,  the  bow  and  stem  both  experience  lift  force  opposed  to  v  and  w,  therefore  Yv  and 

Zw  are  always  negative.  The  terms  Y.  and  are  always  negative  and  have  a  magnitude 
approximately  equal  to  the  displacement  of  the  vessel.  Nv  and  Mw  are  usually  negative, 


but  can  become  positive  if  the  rudder  or  stem  planes  are  very  large.  N.  and  M.  are 
usually  relatively  small  quantities  of  uncertain  sign.15  The  results  shown  in  Table  5  and 
Table  6  match  the  expected  values. 

The  results  contained  in  Table  5  and  Table  6  are  graphically  presented  in  Figure  12 
through  Figure  27.  The  results  have  been  made  non-dimensional  according  to  the  method 
described  in  section  4. 
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Figure  12.  Y ■'  and  Zj  vs  L/Subm  at  Froude  Number  =  0.128 
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Figure  13.  Yv'  and  Zw'  vs  L/Subm  at  Froude  Number  =  0.128 


Added  Mass  vs  L/Subm  at  Froude  Number  =  0.383 


Figure  14.  Y'  and  Zj  vs  L/Subm  at  Froude  Number  =  0.383 


Damping  vs  L/Subm  at  Froude  Number  =  0.383 
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Figure  15.  Yv'  and  Zj  vs  L/Subm  at  Froude  Number  =  0383 


Added  Mass  vs  Froude  Number  at  L/Subm  =  1.277  for 
Sway  and  1.42  for  Heave 
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Figure  16.  Y+  and  Zj  vs  Froude  Number  at  L/Subm  =  1.277  for  Sway  and  1.42  for  Heave 


Damping  vs  Froude  Number  at  L/Subm  =  1.277  for  Sway 

and  1.42  for  Heave 


Figure  17.  Yj  and  Zj  vs  Froude  Number  at  L/Subm  =  1.277  for  Sway  and  1.42  for  Heave 


Added  Mass  vs  Froude  Number  at  L/Subm  =  2.749  for 
Sway  and  2.547  for  Heave 
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Figure  18.  Y-  and  ZM 


vs  Froude  Number  for  L/Subm  =  2.749  for  Sway  and  2.547  for  Heave 


Damping  vs  Froude  Number  at  L/Subm  =  2.749  for  Sway 
and  2.547  for  Heave 


Figure  19.  Yv'  and  Zj  vs  Froude  Number  for  L/Subm  =  2.749  for  Sway  and  2.547  for  Heave 


Figure  12  through  Figure  15  show  that  the  effect  of  submergence  on  the  direct  added 
mass  and  damping  terms  for  sway  and  heave  motion  is  nearly  negligible.  Figure  1 6 
through  Figure  19  show  that  there  is  a  significant  effect  of  speed  on  these  terms.  These 
figures  also  show  that  the  direct  terms  in  sway  and  heave  behave  very  similarly.  The 
cross-term  added  mass  and  damping  coefficients  shown  in  Figure  20  through  Figure  27 
have  different  patterns  of  behavior 


Figure  21.  Nj  and  Mj  vs  L/Subm  for  Froude  Number  =  0.128 
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Figure  22.  AL  and  M w  vs  L/Subm  for  Froude  Number  =  0.383 
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Figure  23.  N'  and  M  '  vs  L/Subm  for  Froude  Number  =  0.383 
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Figure  26.  N x'  and  vs  Froude  Number  for  L/Subm  =  2.749  for  Sway  and  2.574  for  Heave 


The  cross-terms  due  to  sway  and  heave  motion  are  expected  to  be  very  similar  to  each 
other  and  have  similar  dependencies,  however  Figure  20  through  Figure  27  show  that  this 
may  not  be  the  case.  Examination  of  the  measured  moments  listed  in  Table  5  and  Table  6 
show  that  the  measured  yaw  moments  due  to  sway  force  range  from  0.085  to  0.352  N-m 
and  the  measured  pitch  moments  due  to  heave  force  range  from  0.293  to  0.966  N-m.  The 


transducer  is  rated  for  up  to  28  N-m  in  yaw  and  pitch.  This  means  that  the  maximum 
measured  moments  represent  less  than  4%  of  the  range  of  the  transducer.  The  minimum 
moment  represents  only  0.3%  of  the  range  of  the  transducer.  The  abnormal  behavior  of 
the  moment  coefficients  may  be  due  to  insufficient  moment  being  present  to  properly 
deflect  the  transducer. 

The  effects  of  submergence  and  speed  were  determined  numerically  by  performing  a 
least  squares  fit  to  the  data  using  linear  regression.  This  procedure  is  described  in  section 
5.3.4.  The  resulting  equations  for  the  coefficients  are: 

Y '  =  -0.081458  +  0.001774 — - +  0.121754Fr 

Submergence 

Y'  =  -0.016345  +  0.000061 — LenSth - 0.007224Fr 

Submergence 

Z'  =-0.086685  +  0.00091 — - +  0.1 4075 \Fr 

"  Submergence 

Z/  =  -0.013633-0.002679 — length - 0.005671Fr 

Submergence 

M'  =  0.023224  -  0.002943 — - 0.02235 9Fr 

Submergence 

M'  =  0.002825  +  0.000594 — LenSth - o.006407Fr 

Submergence 

N  '  =  -0.00207  +  0.000093 — - 0.016489Fr 

v  Submergence 

N-  =  -0.00089  +  0.000037 — Le-n-^--- - +  0.000172Fr 

Submergence  (36J 

The  quality  of  the  fit  for  these  equations  is  expressed  in  terms  of  the  root-mean-square 

value  of  the  percent  difference  between  the  predicted  and  the  empirical  coefficients. 

Table  7  contains  the  values  for  the  quality  of  fit. 


Table  7.  Quality  of  Fit  for  Heave  and  Sway  Motion  Coefficients 


Coefficient 

Fit 

15% 

y; 

4% 

zj 

17% 

z ; 

9% 

m; 

19% 

m ; 

27% 

< 

9% 

a^; 

9% 

5.4.5.  Forces  and  Moments  Due  to  Body  Angle 

The  forces  and  moments  due  to  body  angle  were  determined  using  data  from 
experiments  performed  on  the  small  scale  model  at  MIT.  Appendix  L  contains  the  list  of 
experiments  performed. 

The  equations  for  the  coefficients  due  to  body  angle  are16 

Y 


X.  =■ 


~^pU2D2a 


Z . =■ 


M.„„  =  ■ 


| pU2D2a 
M 

~ pU2D2La 


(37) 


where 


N 


N 


~pU2DzLa 


D  is  the  diameter  of  the  vehicle 
L  is  the  length  of  the  vehicle 
a  is  the  angle  of  attack. 


The  results  of  these  experiments  are  included  as  Appendix  M.  The  effect  of  the  angles 
is  assumed  to  be  linear  for  small  angles,  so  the  coefficients  were  calculated  for  each 
experiment  and  then  averaged  for  each  combination  of  speed  and  submergence.  Pitch 
and  Yaw  angle  coefficients  were  averaged  separately.  Any  data  point  with  a  calculated 
coefficient  more  than  1.15  standard  deviations  from  the  mean  was  removed  from 
consideration  as  unreliable  data.  1.15  standard  deviations  was  chosen  as  the 
discrimination  point  in  order  to  remove  the  clearly  bad  data  while  retaining  as  much  of 
the  possibly  good  data  as  possible.  The  mean  coefficients  for  each  combination  of 
submergence  and  velocity  are  shown  in  Table  8.  Figure  28  through  Figure  31  illustrate 
the  dependency  of  the  restoring  forces  on  submergence  and  speed. 

It  is  very  important  to  note  at  this  point  that  the  magnitudes  of  the  yaw  and  pitch 
moments  measured  during  these  tests  are  very  small,  on  the  order  of  less  than  1%  of  the 
capacity  of  the  load  cell.  The  data  is  analyzed  and  presented  here,  but  further  work  is 
required  to  determine  the  accuracy  of  the  equations  with  a  more  appropriate  transducer. 
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Table  8.  Mean  Coefficients  for  Force  and  Moment  Due  to  Body  Angle  at  Various  Submergences  and 
Velocities 


Nuv  | 

Submergence 

Velocity 

Nuv 

m 
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0.252 

0.75 

-1.02061 

0.398 

0.75 

-0.99671 

0.543 

0.75 
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1 

-0.80028 

0.398 

1 

-0.73773 

0.543 

1 

-0.58644 

Yuv 

Submergence 

Velocity 

Yuv 

m 

m/s 

0.252 

0.75 

-5.36887 

0.398 

0.75 

-4.85345 

0.543 

0.75 

-2.32929 

0.252 

1 

-3.51475 

0.398 

1 

-2.8487 

0.543 

1 

-1.61492 

Zuw 


Submergence 

Velocity 

Zuw 

m 

m/s 

0.252 

0.75 

1.023955 

0.398 

0.75 

2.269136 

0.543 

0.75 

1.141877 

i  0.252 

1 

1.220249 

0.398 

1 

2.612357 

0.543 

1 

1.801033 

Muw 


Submergence 

Velocity 

Muw 

m 

m/s 

0.252 

0.75 

-0.626 

0.398 

0.75 

-0.59447 

0.543 

0.75 

-0.71149 

0.252 

1 

-0.4012 

0.398 

1 

-0.61916 

0.543 

1 

-0.59916 
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Figure  29.  Zuw  as  a  Function  of  Submergence  for  Two  Speeds 


L/Subm 


+ - ’Fr=  0.290  — ■ — ’Fr=0.386 


Figure  30.  Nuv  as  a  Function  of  Submergence  for  Two  Speeds 


Figure  31.  Muw  as  a  Function  of  Submergence  for  Two  Speeds 

The  data  from  the  restoring  force  and  moment  experiments  was  put  into  a  linear 
regression  model  to  determine  the  dependency  of  those  forces  and  moments  on 
submergence  and  speed.  The  resulting  equations  are 

Yv  =3.035943-11.232476 — —  +  2.399695f — — 1  +15.79518 

Submergence  ^  Submergence  j 


+  1 5.795 188Fr 


Zm.  =  -7.723764  +  9.156976  - 2.365368 

Submergence 

Nuv  =  0.02001 7-1 .452965 — - +  0.317964 

Submergence 


'  Length 
v  Submergence 

Length  ^ 
Submergence  j 


+  4.182974Fr 


+  1.987686Fr 


Muw  =  -1 .1 68727  +  0.127681  — - 0.007568 

Submergence 


Length 

Submergence 


+  1.079278Fr 
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5.4.6.  Sway  Force  and  Yaw  Moment  Due  to  Yaw  Motion 

Experiments  were  performed  with  several  combinations  of  submergence,  velocity,  and 
frequency  of  oscillation  to  determine  the  effect  of  yaw  motion  on  sway  force  and  yaw 
moment.  The  experiments  performed  are  listed  in  Appendix  B.  The  method  of  analysis 
was  similar  to  that  used  for  the  sway  and  heave  motion  tests,  but  in  this  case  there  is 
another  term  due  to  the  angle  of  the  body  as  it  moves  forward  and  oscillates  in  yaw. 

The  equations  are 


and 


Y  =  Yrr  +  Y-r  +  Yuv{-\  pU2D2  a) 

N  =  Nrr  +  N-r  +  Nuv{-\  pU2D2  La) 

Im  (7) 


(39) 


T  =■ 


Y  = 


N  = 


N>  = 


coA 

Re(T)-Tw(~|/C7l72D2a) 
-co2  A 

Im(AQ 

coA 

-co1  A 


(40) 


The  restoring  force  is  a  real  force  and  must  be  subtracted  from  the  measured  real 
hydrodynamic  force  to  calculate  the  force  due  to  the  yaw  motion. 

The  test  conditions  and  results  are  presented  in  Table  10. 

Table  10  Test  Conditions  and  Results  for  Sway  Force  and  Yaw  Moment  Due  to  Yaw  Motion 
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MEM 
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HGSSSIH 

Hz 
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N 

kg 
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KB 

i  □ 

0.252 

IcEEl 
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10 
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B.EWSlrJ 
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MEEMM 
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10 

1.114 
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K>mm 

EBB 
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■UKSKM 

MiZEEM 

10 

■HIM 

-21.2 

WKEES3M 

UliUlfj 

giTiCTZffl 

gifiTi»/;! 

0.543 

ftlcKk— 

1.196 

10 

0.7008 

-0.0670 

gijiffliEii 

ElKUlhUKI 

0.398 

beem 

■IZEM 

10 

WEEESM 

-22.7 

0.3229 

bseeem 

BiIkMiVJ 

EEM1 

0.543 

mem 

bem 

io 

■l  EMM 

-27.8 

■EEEMH 

iTgiTVll 

EEIE3 

graszvi 

5.4.7.  Heave  Force  and  Pitch  Moment  Due  to  Pitch  Motion 

The  analysis  for  heave  force  and  pitch  moment  due  to  pitch  motion  closely  mirrors  the 
analysis  for  sway  force  and  yaw  moment  due  to  yaw  motion.  For  pitch  motion,  the 
equations  are 


Z  =  Zqq  +  Zjfl  +  Zuw(—\pU2D2a) 

M  =  M  q  +  M-q  +  Muw{-\ pU2D2La) 


(41) 


and 
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Im  (Z) 
q  coA 

z 

*"  ^  (42' 

M  =MM1 

q  coA 

Re(M)-MMH,(-j^2£>2Ia) 
q  -a>2A 

The  test  conditions  and  results  are  presented  in  Table  1 1 .  During  these  tests  special 
consideration  was  given  to  removing  the  effects  of  surge  motion  caused  by  the  long  pitch 
arm  of  the  test  apparatus.  In  order  to  remove  these  effects,  the  test  apparatus  was 
oscillated  in  surge  at  the  same  time  pitch  oscillations  occurred.  The  surge  oscillations 
were  180°  out  of  phase  with  the  pitch  oscillations  and  of  such  a  magnitude  as  to  cancel 
the  surge  due  to  pitch. 

Table  11.  Test  Conditions  and  Results  for  Heave  Force  and  Pitch  Moment  due  to  Pitch  Motion 


Pitch  Motion 

Hydrodynamic  Heave  Force,  Z  1 

Submerqence 

Re(Z) 

msssm 

Z9_ 

m 

■rnn 

g»?5iff5!53i 

MK9H 

N 

N 

■IHB 

■ 

■K3H 

■SB 

■  ■ 

0.252 

1.19629 

9 

■cKHcM 

BH'f 

EgHSBI 

0.252 

WKESSM 

WKmm 

9 

-0.4326 

BHI 

ii 

0.543 

■*mi 

aai 

9 

Ki££l!l!l 

KEI2EK 

KEEH9I 

eeebsi 

EEsg 

0.543 

mmmm 

9 

EM 

WKE5HM 

g*i«w[ii 

0.252 

»mnif 

wmm 

9 

■>T£M 

mn TiiiTii 

0.1579 

■»eee» 

gtltlil-VJ 

Pitch  Motion 

Hydrodynamic  Pitch  Moment,  M  | 

Submerqence 

GSIEE31 

gaiiWi 

KIEUJK 

ESSJflZdSj 

■M3M 

WEBm 

ill 

m 

Km 

Hz 

■  ■I3*H;Ekg 

■2SSH 

ISEfflEfSl 

N-m 

N 

_ ks _ _ 

N-m 

wk&uM 

0.252 

■<***■ 

KKEZB 

9 

BBS! 

0.5470 

0.5749 

jgtKiffliEM 

0.252 

BCEUSH 

9 

iKX!Z£9i 

WtEEEBSM 

0.543 

■oftsag 

1.196 

9 

■iWM 

EMI 

■EE22B 

gtity/tyg 

kh 

EEES3 

gifcUfefcEl 

EtltlihidJ 

0.543 

keesb 

9 

BT'Tg 

EMI 

KE23B 

Em 

BiMEHl 

0.252 

Kgsaa 

9 

HMM1 

0.1062 

KKBsBM 

bsb 

»TUlX!« 

eeb 

5.4.8.  Discussion  of  Yaw  and  Pitch  Motion  Results 

The  results  of  yaw  and  pitch  motion  testing  are  graphically  illustrated  in  Figure  32 
through  Figure  47.  Both  submergence  and  speed  were  found  to  be  significant  factors. 
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Added  Mass  vs  L/Subm  at  Froude  Number  =  0.128 
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Figure  32.  N-'  and  M (  '  vs  L/Subm  at  Froude  Number  =  0.128 


Damping  vs  L/Subm.  at  Froude  Number  =  0.128 


0.00000 

-0.00100. 

-0.00200 

-0.00300 

-0.00400 

-0.00500 

-0.00600 

-0.00700 

-0.00800 

-0.00900 

-0.01000 


000- 


-+300- 


2.000- 


2.500 


v  =  -0.001 8x-  0.0035 


y  =  -0.0007x-~0: 


3300 


+  Nr' 
x  Mq' 

- Linear  (Nr') 

- -  Linear  (Mg’) 


Length/Submergence 


Figure  33.  N  '  and  Mg'  vs  L/Subm  at  Froude  Number  =  0.128 


Added  Mass  vs  L/Subm  at  Froude  Number  =  0.383 
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Figure  34.  N-'  and  M(,'  vs  L/Subm  at  Froude  Number  =  0.383 


Damping  vs  L/Subm.  at  Froude  Number  =  0.383 
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Figure 35.  Nr'  and  M vs  L/Subm  at  Froude  Number  =  0.383 
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Figure  36.  N-  and  M ^  vs  Froude  Number  at  L/Subm=  1.277 
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Figure  37.  Nr  and  M  vs  Froude  Number  at  L/Subm=  1.277 


Added  Mass  vs  Froude  Number  for  L/Subm  =  2.749 
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Figure  38.  N /  and  M  /  vs  Froude  Number  at  L/Subm=  2.749 


Damping  vs  Froude  Number  for  L/Subm  =  2.749 
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Figure  39.  Nr  and  Mq  vs  Froude  Number  at  L/Subm=  2.749 


Figure  41.  Yr  and  Zq  vs  L/Subm  at  Froude  Number  =  0.128 
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Figure  43.  Yr  and  Z'q  vs  L/Subm  at  Froude  Number  =  0.383 
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Figure  46.  Y-  and  Z-  vs  Froude  Number  at  L/Subm=  2.749 
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Figure  47.  Yr  and  Z  '  vs  Froude  Number  at  L/Subm=  2.749 


The  data  from  pitch  and  yaw  motion  was  also  analyzed  by  linear  regression  to  yield 
equations  for  the  coefficients.  Those  equations  are 


♦ 


Yr '  =  0.020906  +  0.000403  Len'<th - 0 .0 1 7979Fr 

r  Submergence 

Y\=  0.002324  -  0.000087  — hF!13th - 0.000632Fr 

Submergence 

Z  ’  =  -0.070199  +0.003498 — - +  0.133729Fr 

9  Submergence 

Z  '  =  -0.002666  -  0.000514 — LenSth - +  o.020222Fr 

9  Submergence 

M  '  =  -0.0105 1 5  -  0.000701 — - +  0.030777 Fr 

9  Submergence 

M.  ’  =  0.00014  +  0.000106 — - 0.009854Fr 

9  Submergence 

Nr '  =  -0.006013  -  0.001 173 — —  +  0.012297Fr 

Submergence 

N,  ’  =  -0.000786  -  0.00001 9  — ^  _  0.000582Fr 

Submergence 

Table  12  contains  the  values  for  the  quality  of  fit. 

Table  12.  Quality  of  Fit  for  Pitch  and  Yaw  Motion  Coefficients 


Coefficient 

Quality  of  Fit 

Yr' 

7% 

y; 

46% 

z; 

0% 

0% 

< 

0% 

< 

0% 

K' 

23% 

n; 

8% 

(43) 


6.  Computational  Analysis 

Numerical  method  results  to  accompany  experimental  results  are  very  important.  To 
that  end  the  validated  free  surface  linear  code  for  surface  ships,  SWAN,  was  revised  by 
one  of  its  developers  to  include  submerged  objects.  Unfortunately,  the  results  from  the 
submerged  vehicle  version  of  SWAN  were  too  unreasonable  to  be  either  valid  or 
included  here.  For  example,  in  some  conditions,  the  added  mass  was  predicted  to  be  16 
times  the  actual  displaced  mass.  That  cannot  be  correct.  There  was  not  enough  time  to 
properly  revise,  test,  and  validate  the  new  version  of  this  numerical  method.  Zero 


58 


+ 


# 


forward  speed  evaluations  of  added  mass  and  damping  were  performed  with  the  validated 
code  WAM1T  and  those  results  were  reasonable.  However,  the  interest  for  this  research 
is  for  an  underwater  vehicle  with  forward  speed. 

7.  Conclusion 

The  analyses  performed  for  this  research  have  drawn  upon  a  limited  number  of 
experiments  with  a  great  deal  of  uncertainty  in  the  quality  of  data.  The  limited  number  of 
experiments  is  a  result  of  very  significant  delays  experienced  in  designing  and  testing  the 
test  apparatus.  Although  limited  data  exists  that  is  known  to  be  valid,  a  very  large 
amount  of  data  exists  from  earlier  experiments.  These  earlier  experiments  may  have  had 
interaction  between  the  transducer  and  the  vehicle  shell.  Although  the  geometry  of  the 
model  indicates  no  interference,  there  may  have  been  interference  due  to  shell  distortion 
and  corrosion  products.  Also,  some  of  the  uncertainty  in  the  quality  of  early  data  results 
from  the  measurement  of  very  small  forces  and  moments  that  represent  a  very  small 
fraction  of  the  range  of  the  transducer.  The  method  of  research  and  analysis  is  sound  and, 
at  the  very  least,  provides  a  roadmap  for  conducting  the  research  in  the  future. 

The  results  associated  with  the  sway  force  due  sway  motion  and  heave  force  due  to 
heave  motion  are  very  similar  and  indicate  that  the  quality  of  that  data  is  very  good.  The 
effects  of  variation  in  submergence  and  speed  on  those  forces  can  also  be  expected  to  be 
good  results. 

The  proximity  to  the  free  surface  has  relatively  little  effect  on  forces  due  to  motions  or 
control  surface  deflections.  This  indicates  that  with  minor  control  system  alterations, 
operation  in  shallow  water  is  probably  feasible.  Of  course  further  work  on  proximity  to 
the  bottom  is  necessary. 

The  causes  for  the  variation  in  the  results  with  Froude  number  and  submergence  are 
not  well  understood.  There  are  various  theories  that  discuss  the  interaction  of  the  flow 
stream  with  the  free  surface  and  the  bottom  that  may  be  useful,  but  true  understanding  of 
the  results  of  these  analyses  must  be  saved  for  future  work. 

8.  Future  Work 

In  this  thesis  we  have  found  the  effects  of  shallow  water  and  submergence  on  many  of 
the  dominant  terms  in  the  maneuvering  equations.  The  reason  for  these  variations  is  also 
yet  to  be  determined.  In  the  very  near  future  the  data  from  the  experiments  conducted 
early  in  this  research  and  not  known  to  be  reliable  will  be  analyzed.  If  that  data  is  found 
to  be  reliable,  the  results  reported  here  will  be  expanded  and  modified  to  include  that  data 
and  provide  a  better  description  of  the  effects  of  submergence  and  speed. 

This  work  indicates  that  the  proximity  to  the  free  surface  has  relatively  little  effect  on 
forces  due  to  motions  or  control  surface  deflections.  This  indicates  that  with  minor 
control  system  alterations,  operation  in  shallow  water  is  probably  feasible.  Of  course 
further  work  on  proximity  to  the  bottom  is  necessary.  Much  of  this  is  available  from  the 
work  of  William  Ramsey.17 

Additional  future  work  will  involve  finding  a  numerical  method  that  can  accurately 
predict  the  phenomena  studied  here.  The  results  of  this  and  similar  research  will  be  used 
to  validate  that  method  for  widespread  use  in  predicting  this  kind  of  behavior. 
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Appendix  A.  Full  Scale  Experiments  Performed  at  the 
United  States  Naval  Academy  to  Determine  the  Effects  of  Body 


Angle  and  Control  Surface  Deflection 

NAVAL  ACADEMY  REMUS  NO  WAVE  STEADY  FORCE  TEST  RUNS 
_ (depth  is  measured  to  top  of  vehicle  at  strut  center) 
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• 

• 
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Rudder 

Angle 

(deg) 
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Plane 

Angle 

(deg) 

Pitch 

Angle 

(deg) 

Yaw 

Angle 

(deg) 

Speed 

(ft/s) 

Depth 

(meters) 

V/ 

sqrt(gL) 

depth/L 

N1 

0 

0 

0 

0 

6.8 

0.762 

0.526 

0.481 

N2 

0 

0 

0 

8 

3.4 

0.762 

0.263 

0.481 

N3 

0 

0 

8 

0 

6.8 

0.762 

0.526 

0.481 

N4 

0 

7 

0 

0 

6.8 

0.762 

0.526 

0.481 

N5 

7 

0 

0 

0 

3.4 

0.762 

0.263 

0.481 

N6 

0 

0 

0 

0 

1.69 

0.381 

0.131 

0.240 

N7 

0 

0 

0 

0 

HI 

0.381 

0.263 

0.240 

N8 

0 

0 

0 

0 

6.8 

0.381 
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0.240 

N9 

7 

0 

0 

0 

6.8 
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0.240 

N10 

0 

7 

0 

0 

6.8 
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0.526 
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Nil 
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0.381 
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N20 

0 

7 

0 

0 

6.8 
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0.120 

N21 
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N22 

7 

0 

0 

0 

6.8 

0.1905 

0.526 

0.120 

N23 

7 

0 

-8 

0 

6.8 

0.1905 

0.526 

0.120 

N24 

7 

7 

0 

0 

6.8 

0.1905 

0.526 

0.120 
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Appendix  B.  Model  Scale  Experiments  Performed  at  MIT  to 
Determine  the  Forces  and  Moments  Due  to  Unsteady  Motion 


Velocity 
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0.333 
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lOSClOOl  0.795  I  0.252  | _ |  I  | _ | _ | _ |  10  j  0.4028  |  0.333  | 


Appendix  C.  AutoanalyzeXls.m 

%  AutoanalyzeXls.m 
%  Erik  Oiler,  2003 

%  This  program  uses  other  programs  to  automatically,  analyze  all  the  data 
%  files. 

%  Operational  Overview 

%  1 .  User  starts  this  program  to  perform  analysis  of  experimental  data 
%  files. 

%  2.  User  selects  the  series  of  experiments  to  be  analyzed. 

%  3.  The  program  opens  an  excel  file  called  "MIT  Test  Plan.xls"  and 
%  imports  the  data  from  the  worksheet  for  the  selected  series.  The 
%  matrix  seriesnum  contain  the  numerical  data  from  the  worksheet 
%  and  the  matrix  seriestext  contains  the  text  data  from  the 
%  worksheet. 

%  4.  The  program  runs  the  appropriate  analysis  program  for  each  test 
%  series.  For  steady  force  tests,  the  analysis  program  is 
%  "AnalyzemodXlsSF.m".  For  all  other  tests,  the  analysis  program 
%  is  "AnalyzemodXls.m". 

%  5.  The  called  analysis  program  analyzes  the  raw  data  and  writes  the 
%  will  display  that  the  data  file  does  not  exist. 

%  6.  When  all  data  files  listed  in  the  test  plan  have  been  analyzed, 

%  file. 

%  7.  The  following  files  must  be  in  the  same  directory: 

%  "AutoanalyzeXls.m" 

%  "AnalyzemodXlsSF.m" 

%  "AnalyzemodXls.m" 

%  "MIT  Test  Plan.xls" 

%  All  Data  Files  to  be  analyzed. 

%  Initialize  the  workspace  by  clearing  all  variables  and  closing  all 
%  windows, 
close  all; 
clear  all; 


%  Determine  which  set  of  tests  to  analyze. 

fprintf('  1 :  Horizontal  Plane  \n') 

fprintf('  2:  Vertical  Plane  \n') 

fprintf('  3:  Pure  Sway  \n') 

fprintf('  4:  Pure  Heave  \n') 

fprintf('  5:  Pure  Pitch  \n') 

frrintff  6:  Pure  Yaw  \n') 

fprintf('  7:  Mass  Matrix  \n') 


fprintf('  8:  Steady  Force  \n') 
fprintf('  9:  Inertial  Calculation  Checks  \n') 
fprintf('10:  Miscellaneous  \n') 
fprintf('l  1 :  Oscillation  Tests  \n') 

series  =  input('Which  Test  Series?  \n'); 

%  Select  the  worksheet  in  "MIT  Test  Plan.xls"  and  the  output  file  based 
%  upon  the  test  series, 
switch  series 
case  1 

worksheet  =  'Hor  Plane'; 
outputfile  =  'HorPlaneOut.txt'; 
case  2 

worksheet  =  'Compensating  Vert  Plane'; 
outputfile  =  'VertPlaneOut.txt'; 
case  3 

worksheet  =  'Pure  Sway'; 
outputfile  =  'SwayOut.txt'; 
case  4 

worksheet  =  'Pure  Heave'; 
outputfile  =  'HeaveOut.txt'; 
case  5 

worksheet  =  'Pure  Pitch'; 
outputfile  =  'PitchOut.txt' ; 
case  6 

worksheet  =  'Pure  Yaw'; 
outputfile  =  'Y  awOut.txt'; 
case  7 

worksheet  =  'Mass  Matrix  Tests'; 
outputfile  =  'MassTestsOuUxt'; 
case  8 

worksheet  =  'Steady  Force'; 
outputfile  =  'SteadyForceOut.txt'; 
case  9 

worksheet  =  'Mass  Matrix  Tests'; 
outputfile  =  'InertiaCalcTestsOut.txt'; 
case  10 

worksheet  =  'Miscellaneous'; 
outputfile  =  'MiscellaneousOut.txt'; 
case  1 1 

worksheet  =  'Oscillation'; 
outputfile  =  'OscillationOut.txt'; 
end 

%  Import  the  test  filenames  and  test  conditions  from  the  selected 


%  worksheet. 

[seriesnum,seriestext]  =  xlsread('MIT  Test  Plan', worksheet); 
numtests  =  size(seriesnum,l)-2; 

if  series  =1 1  %  For  an  unknown  reason,  the  oscillation  test  worksheet 
%  imports  differently  from  the  other  worksheets, 
numtests  =size(seriesnum,  1 ); 
for  testindex  =  1  rnumtests 

run(testindex)  =  (seriestext(testindex+2,3)); 

DandS(testindex,l)  =  seriesnum(testindex,l); 

DandS(testindex,2)  =  seriesnum(testindex,2); 
end 
else 

for  testindex  =  1  rnumtests 

run(testindex)  =  (seriestext(testindex+2,3)); 

DandS(testindex,l)  =  seriesnum(testindex+2,l); 
DandS(testindex,2)  =  seriesnum(testindex+2,2); 
end 
end 

%  Display  the  first  and  last  files  to  be  analyzed. 
firsttest=char(run(  1 )) ; 
lasttest=char(run(numtests)) ; 

fprintf('Autoanalyze  will  proceed  from  %s  to  %s  \n',  firsttest,lasttest) 

%  Open  the  output  file, 
warning  off 
delete(outputfile) 
warning  on 

manyrowsfid=fopen(outputfile,'a'); 

%  Display  which  file  is  being  processed  and  process  the  test  file.  Send 
%  the  test  condition  data  that  can  not  be  extracted  for  the  data  file, 
for  fileindex  =  1  rnumtests 
fname  =  char(run(fileindex)); 
fprintf('Processing  %s\n', fname) 

Depth  =  DandS(fileindex,l); 

Submergence  =  DandS(fileindex,2); 
if  series  ==  8  %  Steady  force  tests. 

Yaw  Angle  =  seriesnum(fileindex+2,3); 

PitchAngle  =  seriesnum(fileindex+2,4); 

Rudder  =  0; 

StemPlanes  =  0; 
analyzemodxlssf; 
else  %  All  unsteady  motion  tests. 


Y  aw  Angle  =  0; 
PitchAngle  =  0; 
Rudder  =  0; 
StemPlanes  =  0; 
analyzemodxls 
end 


end 

%  Display  that  processing  is  complete  and  close  the  output  file. 
fprintf('Processing  Complete'); 
status  =  fclose(manyrowsfid); 


Appendix  D.  Sample  Data  File 

This  is  the  first  portion  of  the  data  from  experiment  OSC07.  OSC07  involved 
oscillations  in  sway  with  forward  velocity.  The  two  rows  of  text  at  the  top  were  added  by 
the  author  for  clarification  and  are  not  actually  part  of  the  data  file.  Not  shown  are  the 
additional  rows  at  the  end  of  the  data  file  which  contain  the  ordered  parameters  for  the 


experiment. 

|  Forces  and  Moments  j 

Wave  Height 
Channels  (Not  Used) 

j  Position  Data  in  Controller  Counts  j 

i  Time  Data  I 

Alternate  Time  Data 

mmmm 

|X  |Y 

|Roll  !  Pitch  i 

|Time  Data|Da1 

te  Data|DeltaT 

Time  Date|DettaT 

0.817 

-0.031 

0.086 

0.086 

0.102 

-0.115 

0.002 

0.001 

165 

316744 

798231 

503551 

13916 

37142.76 

37725 

0 

0.027 

0 

0.322 

-0.099 

0.005 

0.021 

-0.456 

-0.026 

0 

0.001 

166 

319080 

798231 

503551 

13916 

37142.8 

37725 

0.039063 

0,073 

0.046 

0.667 

-0.057 

0.005 

0.072 

-0.796 

0.002 

-0.001 

0.001 

227 

321424 

798231 

503551 

13916 

37142.84 

37725 

0.039063 

0.116 

0.043 

0.361 

-0.006 

0.014 

-0.091 

0.784 

-0.079 

0.002 

0.001 

416 

323714 

798231 

503551 

13916 

37142.88 

37725 

0.042969 

0.156 

0.04 

-0.052 

-0.016 

0.005 

0.149 

-0492 

-0.045 

0.001 

0.002 

836 

326137 

798230 

503551 

13916 

37142.92 

37725 

0.039063 

0.198 

0.042 

0.413 

-0.036 

0 

0.584 

-0  238 

0.012 

0.001 

0.003 

1350 

327779 

798231 

503551 

.13916 

37142.96 

37725 

0.039063 

0,241 

0.043 

0.638 

-0.074 

0.003 

0.347 

-0.371 

-0.079 

0.001 

0.001 

2548 

330343 

798230 

503551 

13916 

37143 

37725 

0.039063 

0.269 

0028 

0.294 

-0.086 

-0.005 

0.489 

0.129 

-0.111 

0 

0.001 

4159 

332879 

798232 

503551 

13916 

37143.0 4 

37725 

0.039063 

0.31 

0.041 

-0.027 

-0.046 

0.005 

0.295 

-0.067 

-0.049 

0.001 

0 

6398 

335523 

798230 

503551 

13916 

37143.08 

37725 

0.042969 

0.352 

0.042 

0.229 

-0.101 

0.008 

-1.175 

0.169 

-0.065 

0.001 

0.001 

9262 

338054 

798231 

503551 

13916 

37143.12 

37725 

0.039063 

0.395 

0.043 

-0.008 

0.011 

0.013 

-1.092 

0.119 

-0.004 

0 

0 

12747 

340486 

798230 

503551 

13916 

37143.16 

37725 

0.039063 

0.437 

0.042 

0.286 

0.024 

0.004 

0.519 

-0.023 

0.033 

0 

0.001 

16652 

342728 

798231 

503551 

13916 

37143.2 

37725 

0.039063 

0.477 

0.04 

0.345 

-0059 

-0.005 

0.794 

-0.127 

-0.118 

-0  006 

0 

21057 

344921 

798230 

503551 

13916 

37143.24 

37725 

0.042969 

0.518 

0.041 

0.043 

-0.167 

0.004 

0.646 

0.224 

-0  166 

0 

0.001 

25744 

347005 

798231 

503551 

13916 

37143.28 

37725 

0.039063 

0.558 

0.04 

-0.057 

0.042 

0.011 

0.361 

-0.055 

0.07 

0 

0.002 

29136 

348342 

798231 

503551 

13916 

37143.32 

37725 

0.039063 

0.587 

0.029 

-0.341 

-0.067 

0.015 

-0.454 

0.676 

0.085 

0 

0.002 

34368 

350165 

798230 

503552 

13916 

37143.36 

37725 

0.042969 

0.629 

0.042 

-0.234 

0.013 

0.011 

-0.067 

0.159 

0026 

0.001 

0.007 

39725 

351745 

798232 

503551 

13916 

37143.4 

37725 

0.039063 

0.672 

0.043 

-0.069 

0.1  B4 

0.015 

0.009 

-0.448 

-0.033 

0 

0.002 

45221 

353150 

798231 

503551 

13916 

37143.44 

37725 

0.039063 

0.714 

0.042 

0208 

-0.013 

0.031 

-0.441 

0.257 

-0029 

0.001 

0.001 

50493 

354288 

798231 

503551 

13916 

37143.48 

37725 

0.039063 

0.754 

0.04 

-0.189 

-0.152 

0.012 

0.088 

-0.165 

0.004 

0.001 

0.001 

55941 

355252 

798232 

503551 

13916 

37143.52 

37725 

0.042969 

0.797 

0.043 

-0.843 

-0.053 

0.027 

0.165 

0.B34 

0.084 

-0.001 

0 

61351 

355897 

798232 

503551 

13916 

37143.56 

37725 

0.039063 

0.839 

0042 

-0.409 

-0.003 

0.017 

0.167 

-0.171 

0  021 

0.001 

0 

64685 

356137 

798231 

503551 

13916 

37143.6 

37725 

0.039063 

0.88 

0.041 

-0  109 

0.034 

0.014 

0.562 

0.097 

0.039 

0.001 

0.001 

70086 

356319 

798231 

503551 

13916 

37143.64 

37725 

0.039063 

0.908 

0.028 

-0.173 

-0.03 

0.01 

-0.12 

0.203 

-0.003 

0.001 

0 

75293 

356303 

798231 

503551 

13916 

37143.68 

37725 

0.042969 

0.949 

0  041 

-0.205 

-0.018 

0.007 

-0.403 

-0.553 

0.049 

0 

0 

80797 

356089 

798231 

503551 

13916 

37143.72 

37725 

0.039063 

0.991 

0.042 

-0.698 

-0.089 

0.013 

0.284 

0.863 

0.061 

0.001 

0.001 

86001 

355606 

798232 

503551 

13916 

37143.76 

37725 

0.039063 

1.031 

004 

-0.97 

0.019 

0.015 

0.106 

0.113 

-0.007 

0.001 

0.001 

91428 

354722 

798232 

503551 

13916 

37143.8 

37725 

0.039063 

1.074 

0.043 

-0.219 

-0.052 

0.012 

0.192 

-0.204 

0.035 

0 

0 

96569 

353570 

798231 

503551 

13916 

37143.84 

37725 

0.042969 

1.114 

0.04 

-0.267 

0.033 

0.014 

0.26 

0.544 

0.066 

0 

0.002 

101985 

352184 

798231 

503551 

13916 

37143.88 

37725 

0.039063 

1.157 

0.043 

-0.382 

0.05 

0.065 

-0.104 

-0.105 

-0.024 

0 

0 

107466 

350626 

798231 

503551 

13916 

37143.92 

37725 

0.039063 

1.199 

0.042 

-0.141 

-0.035 

0.012 

-0.241 

-0.123 

0.036 

0.001 

0 

111116 

349500 

798231 

503551 

13916 

37143.96 

37725 

0.039063 

1.228 

0.029 

-0.668 

0.002 

0.019 

0.342 

0.669 

0.02 

0001 

0 

116324 

347693 

796231 

503551 

13916 

37144 

37725 

0.042969 

1.268 

0.04 

-0.688 

0.077 

0.014 

0.321 

-0.271 

0.059 

0 

0.001 

121748 

345560 

798231 

503551 

13916 

37144.04 

37725 

0.039063 

1.311 

0.043 

-0.136 

-0  047 

-0.004 

0.218 

0.261 

-0.003 

0.001 

0.001 

127162 

343249 

798231 

503551 

13916 

37144.08 

37725 

0.039063 

1.353 

0.042 

-0.19 

0.063 

0.003 

0.39 

-0.264 

0.069 

0.001 

0.001 

132580 

340876 

798231 

503551 

13916 

37144.12 

37725 

0.039063 

1.395 

0.042 

-0.062 

-0.066 

0.008 

-0.456 

-0.179 

0.039 

0001 

0.001 

137793 

338568 

798231 

503551 

13916 

37144.16 

37725 

0.042969 

1.436 

0.041 

-0.208 

-0.041 

0.011 

-0.45 

-0.086 

-0018 

0 

0.001 

143271 

336076 

798231 

503551 

13916 

37144.2 

37725 

0.039063 

1.478 

0  042 

-0.6 

-0.094 

0.013 

0.581 

0.243 

0.072 

0.001 

0.001 

147161 

334225 

798232 

503551 

13916 

37144.24 

37725 

0.039063 

1.509 

0.031 

-0.333 

0.107 

0.009 

-0.048 

0  169 

-0.048 

'  0 

0 

152338 

331709 

798232 

503551 

13916 

37144.28 

37725 

0.039063 

1.563 

0.054 

-0.061 

-0154 

0.01 

-0.065 

-0.001 

-0.025 

-0.003 

0.001 

157760 

329060 

798231 

503551 

13916 

37144.32 

37725 

0.042969 

1.592 

0.029 

0.509 

0.079 

0.001 

0.307 

-0  365 

-0.025 

0.002 

0 

163177 

326519 

798231 

503551 

13916 

37144.36 

37725 

0.039063 

1.634 

0.042 

0.11 

-0.038 

0.001 

-0.117 

-032 

0  117 

0 

0.002 

168907 

323971 

798231 

503551 

13916 

37144.4 

37725 

0.039063 

1.679 

0.045 

-0.183 

-0.019 

0.015 

-0.508 

-0  109 

0.009 

0.001 

-0.001 

174105 

321717 

798232 

503551 

13916 

37144.44 

37725 

0.039063 

1.719 

0.04 

-0.264 

-0.197 

0.014 

0.302 

0.189 

-0.049 

0.001 

0.001 

179560 

319427 

798232 

503551 

13916 

37144.48 

37725 

0.042969 

1.761 

0042 

-0.124 

-0.011 

0.008 

0.177 

0.261 

-0.152 

0.001 

0 

183179 

317973 

798232 

503551 

13916 

37144.52 

37725 

0.039063 

1.79 

0.029 

0.584 

-0.259 

0.022 

0.136 

-0.581 

0.012 

0  001 

0 

188340 

316043 

798231 

503551 

13916 

37144.56 

37725 

0.039063 

1.83 

0.04 

0.397 

0.078 

0.005 

0.792 

0.021 

0  092 

0 

0.001 

193773 

314203 

798231 

503551 

13916 

37144.6 

37725 

0.039063 

1.873 

0.043 

0.073 

-0.127 

0.005 

-0.214 

-0.107 

-01 

0 

0.001 

198972 

312694 

798231 

503551 

13916 

37144.64 

37725 

0.042969 

1.913 

004 

0.153 

-0.039 

0.013 

-0.529 

0.065 

-012 

-0  001 

0.002 

204427 

311307 

798232 

503551 

13916 

37144.68 

37725 

0.039063 

1.956 

0.043 

0.091 

-0.033 

0.012 

-0.027 

-0.299 

0073 

0001 

0.001 

209887 

310112 

798232 

503551 

13916 

37144.72 

37725 

0.039063 

1.998 

0.042 

0.572 

-0.008 

0.023 

-1.172 

-0.164 

-0.008 

0.002 

0.001 

215065 

309219 

798232 

503551 

13916 

37144.76 

37725 

0.042969 

2.038 

004 

0.474 

-0.261 

0.019 

-0.152 

0.118 

-0.047 

0.001 

0 

220485 

308554 

798232 

503551 

13916 

37144.8 

37725 

0.039063 

2.081 

0.043 

0.308 

0.067 

0.012 

0.431 

-0.427 

0.049 

0.001 

0 

223845 

308301 

798232 

503551 

13916 

37144.84 

37725 

0.039063 

2.107 

0.026 

0.676 

-0.112 

0.01 

0.179 

-0.08 

-0.019 

0.001 

0 

229289 

308112 

798231 

503551 

13916 

37144.88 

37725 

0.039063 

2.15 

0  043 

0.12 

0.135 

0.017 

0.102 

0.201 

-0069 

0.001 

0 

234754 

308141 

798231 

503551 

13916 

37144.92 

37725 

0.042969 

2.192 

0.042 

0.345 

-0.029 

0.014 

0.504 

-0.113 

-0.03 

0.002 

0.001 

239953 

308375 

798231 

503551 

13916 

37144.96 

37725 

0.039063 

2.233 

0.041 

0.592 

-0.043 

0.011 

-0.115 

-0.098 

-0.062 

0.001 

0.001 

245392 

308866 

798231 

503551 

13916 

37145 

37725 

0.039063 

2.275 

0.042 

0.997 

-0.122 

0.026 

0.062 

-0.601 

-0.046 

-0.001 

0 

250823 

309761 

798231 

503551 

13916 

37145.04 

37725 

0.039063 

2.317 

0.042 

0.852 

0.007 

0.019 

0.266 

0.224 

-0  083 

-0.001 

0.001 

256250 

310980 

798231 

503551 

13916 

37145.08 

37725 

0.039063 

2.36 

0.043 

0.216 

-0.065 

0.013 

-0.168 

-0.435 

0.005 

0.001 

0.002 

261434 

312375 

798231 

503551 

13916 

37145.12 

37725 

0.042969 

2.4 

0.04 

0.158 

-0.005 

0.02 

-0.38 

0.164 

-0.049 

-0.021 

0.003 

265073 

313395 

798231 

503551 

13916 

37145.16 

37725 

0.039063 

2  429 

0029 

0.509 

-0.084 

0.013 

0.315 

-0.457 

-0.018 

-0.001 

0.001 

270532 

315070 

798231 

503551 

13916 

37145.2 

37725 

0.039063 

2483 

0054 

0.697 

-0.043 

0.016 

0.169 

0.255 

-0115 

0.001 

0.003 

275725 

316843 

798232 

503551 

13916 

37145.24 

37725 

0.039063 

2.512 

0.029 

0.476 

-0.116 

0.016 

-0.246 

-0.355 

-0.014 

0 

V  0 

280904 

318874 

798232 

503551 

13916 

37145.28 

37725 

0.042969 

2.552 

0.04 

0.748 

-0.047 

0.021 

-0.064 

-0.639 

0.021 

0.001 

0 

286336 

321203 

798232 

503551 

13916 

37145.32 

37725 

0.039063 

2.594 

0.042 

0.305 

-0.001 

0.051 

-0.252 

0825 

-0  066 

0 

-0.001 

291770 

323613 

798232 

503551 

13916 

37145.36 

37725 

0.039063 

2.637 

0.043 

-0.091 

0.014 

0.012 

-0.507 

-0  449 

•0.05 

0.001 

0.001 

297219 

326056 

798232 

503551 

13916 

37145  4 

37725 

0.042969 

2.679 

0.042 

0.382 

-0.016 

0.021 

0.265 

-0.11 

-0.031 

0.001 

0.001 

300600 

327567 

798232 

503551 

13916 

37145.44 

37725 

0.039063 

2.72 

0.041 

0.648 

-0.076 

0012 

0.397 

-0.193 

-0.112 

0.001 

0 

306050 

330104 

798232 

503551 

13916 

37145.48 

37725 

0.039063 

2.762 

0.042 

0.426 

-0.126 

0.025 

0.019 

-0.047 

-0.095 

0 

0.001 

311488 

332753 

798232 

503551 

13916 

37145.52 

37725 

0.039063 

2.791 

0.029 

■  0.012 

-0.131 

0,016 

0.711 

-0  134 

-0.025 

0 

0.001 

316927 

335402 

798232 

503550 

13916 

37145.56 

37725 

0.042969 

2.833 

0.042 

0.175 

0.014 

0.007 

0.311 

0.11 

-0.062 

0.001 

0.001 

322098 

337845 

798232 

503551 

13916 

37145.6 

37725 

0.039063 

2.874 

0.041 

-0.033 

-0.119 

-0.018 

-0181 

0.614 

•0.076 

0.001 

0.001 

327555 

340282 

798232 

503551 

13916 

37145.64 

37725 

0.039063 

2.916 

0.042 

Appendix  E.  AnalyzemodXlsSF.m 

%  AnalyzemodXlsSF.m 
%  Erik  Oiler,  2003 

%  Performs  analyses  for  Steady  Force  Tests  only 

%  This  program  is  designed  to  be  called  by  AutoAnalyzeXls.m  and  requires 
%no  manual  intervention. 

%  Initialize  matrices. 

[Data]  =  0; 

[Datain]  =0; 

%  Get  the  input  data  and  find  the  length  of  the  file, 
infname  =  strcat(fname,'.xls'); 

existencecheck  =  exist(infname);  %  Determines  if  the  input  file  exists, 
if  existencecheck  =0 

fprintf('%s  does  not  exist.W, infname) 
return 
end 

Datain  =  xlsread(infhame);  %  Reads  the  input  file. 

NumLines  =  size(Datain,l)-4;  %  Gets  the  number  of  data  samples  in 
%  the  input  file. 

%  Gets  the  date  and  converts  from  EXCEL  to  MATLAB  format. 
Date=datestr(Datain(l,15)-36525,2); 

%  Build  a  date-time  string  to  ensure  the  correct  gains  are  applied. 

%  Add  693960  to  the  date  to  get  number  of  days  from  0000. 

%  Divide  the  number  of  minutes  by  86400  to  get  fractional  days. 

%  MATLAB  date  serial  numbers  are  in  the  form: 

%  days  since  0000. fraction  of  a  day 
DateNum=Datain( 1 , 1 5)+693 960+Datain(  1 , 1 4)/86400; 

%  Deterimine  the  ordered  velocity. 

XSpd  =  Datain(NumLines+3,10); 
if  isnan(XSpd)=l 
XSpd  =  0; 
end 

%  Determine  the  ordered  distance  of  travel. 

XDist  =  Datain(NumLines+3,l  1); 
if  isnan(XDist)==T 


XDist  =  0; 
end 

%  Determine  the  ordered  length  of  time  of  the  experiment. 

Time  =  Datain(NumLines+4,5); 

%  Determine  actual  sample  frequency  based  upon  the  measured  interval 
%  between  data  points.  Necessary  because  the  control  software  did  not 
%  always  sample  at  the  ordered  sample  rate. 

SampleFreq  =  l/mean(Datain(l:NumLines,16)); 


%  Drop  the  first  1 .2  seconds  of  data  to  account  for  acceleration.  At 

#  %  the  sample  frequency  of  25  Hz,  drop  the  first  30  points.  Shift  the 
%  other  points  up  in  the  array. 

dropgap  =  round(  1 ,2*SampleFreq); 
for  I  =  1  :NumLines-dropgap; 

Datain(1, 1 : 1 6)=Datain(I+dropgap,  1:16); 

#  end 

NumLines  =  NumLines  -  dropgap; 

%  Drop  data  recorded  after  the  model  stopped  moving  along  the  X-axis 
^  %  for  non-Inertial  tests. 

NumLines2=NumLines; 
for  1  =  NumLines- 1:- 1:75 

if  (Datain(I,9)==Datain(I+l,9))  %Looks  for  constant  X  position. 

NumLines2  =  I; 
end 

®  end 

NumLines  =  NumLines2; 

%  Determine  the  time  of  the  first  data  point.  Used  for  converting 
%  from  time  past  midnight  to  time  of  run. 

#  TimeStart  =  Datain(l,14); 

%  The  inverse  sensitivity  matrix  (B).  The  sensitivity  matrix  used 
%  depends  on  the  when  the  experiment  was  performed.  Ealry  experiments 
%  used  a  different  load  cell  than  later  ones. 


if((datenum(Date)>=datenum(2002,07,17))&... 
(datenum(Date)<datenum(2002, 1 0,3 1 ))) 

[B]  =  [  0.3901  0.0029  -0.0071  -0.0010-0.0024  -0.0016; 
0.0023  0.3887  0.0016  0.0025-0.0039  0.0019; 
0.0139  0.0129  1.5006  -0.0002  -0.0215  -0.0018; 
0.0000-0.0001  0.0016  0.0069  0.0000-0.0001; 
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-0.0001  0.0000  0.0018  0.0000  0.0070  0.0000; 

-0.0003  -0.0002  0.0000  0.0000  0.0000  0.0108]; 
else  %valid  on  and  after  Oct  3 1 ,  2002 

[B]  =  [  0.3856  0.0020  -0.0016  -0.0008  0.0001  -0.0030; 

0.0023  0.3811  -0.0040  0.0012-0.0034  0.0031; 

0.0093  0.0024  1.5109  0.0068-0.0231  -0.0014; 

0.0001  0.0000  0.0013  0.0069  0.0000  0.0000; 

0.0000  0.0000  0.0007  0.0000  0.0069-0.0001; 

-0.0003  0.0001  0.0000  0.0000  0.0000  0.0106]; 
end 

%  Establish  Gains  and  Excitation  Voltage 
if  (DateNum>73 1434.61458)  %  August  6, 2002  1445. 

[Gain]  =  [4000;4000;1000;1000;1000;4000];  %  y  z  x  pitch  yaw  roll 
[Vexc]  =  [10;  10;  2.5;  10;  10;  10]; 
end 

%  Calculate  the  Conversion  Factors  (CF)  between  voltage  and  force 
%  or  moment 

CFtemp  =  Gain.*Vexc*10A-6; 

CF  =  zeros(6,6); 

CF(1,1)=  CFtemp(l); 

CF(2,2)=  CFtemp(2); 

CF(3,3)=  CFtemp(3); 

CF(4,4)=  CFtemp(4); 

CF(5,5)=  CFtemp(5); 

CF(6,6)=  CFtemp(6); 

%  Establish  multipliers  to  convert  from  controller  counts  to  MKS  units. 
%  Multipliers  are  different  for  different  time  intervals  due  to  system 
%  upgrades. 

%  Between  0800  Oct  24, 2002  and  0800  Dec  1,  2002. 
if  ((DateNum>73 1 5 1 3.33333)&(DateNum<73 1551.33333)) 
xfactor  =  3850; 
yfactor  =  2410; 
zfactor  =  2114; 
pitchfactor  =  972; 
yawfactor=  1818; 

elseif  DateNum>73 155 1.33333  %  After  0800  Dec  1,  2002 
xfactor  =  3850; 
yfactor  =  2410; 
zfactor=  2114; 
pitchfactor  =155; 
yawfactor=  1818; 
else  %  Before  Oct  24,  2003 


xfactor  =  3850; 
yfactor  =  2410; 
zfactor  =  4921; 
pitchfactor  =  1111; 
yawfactor  =  1025; 
end 

%  Convert  from  voltages  to  forces 
for  I  =  1  :  NumLines 

Datain(I,l:6)  =  (CFM*B*Datain(I,l:6)')'; 

Datain(I,9)  =  Datain(I,9)  /xfactor; 

Datain(I,10)  =  Datain(I,10)  /  yfactor; 

Datain(I,l  1)  =  Datain(I,l  1)  /  zfactor; 

Datain(I,12)  =  Datain(I,12)  /  yawfactor; 

Datain(I,13)  =  Datain(I,13)  /  pitchfactor; 

Datain(I,14)  =  Datain(I,14)  -  TimeStart; 
end 

%  Find  the  mean  force 
for  J=l:6 

MeanForce(J)  =  real(mean(Datain(l  :NumLines,J))); 
end 

%  PROCESSING  STAGE 

%  Shift  the  origin  of  the  coordinate  system  from  the  origin  of  the  load 
%  cell  to  vessel  amidships, 
shiftlength  =  0.0522;  %  meters 

MeanForce(4)  =  MeanForce(4)  -  MeanForce(2)  *  shiftlength;  %  Pitch 
MeanForce(5)  =  MeanForce(5)  +  MeanForce(l)  *  shiftlength;  %  Yaw 

%  Print  the  results  to  a  common  row  output  file. 
fprintf(manyrowsfid,'%s  \t  %s\t  %6.2f,fname, Date, Depth); 
fprintf(manyrowsfid,'\t%6.2f  \t%5 .  If  \t  %5 .2f\t%5 . 1  f ,. .  . 

Submergence, XSpd,XDist, Time); 
fprintf(manyrowsfid,'  \t  %2.0f  \t%2.0f\t  %3.1f\t  %3.1f,... 

PitchAngle, Yaw  Angle, StemPlanes, Rudder); 
fprintf(manyrowsfid,'  \t  %7.4f,... 

MeanForce(3),MeanForce(l),MeanForce(2),MeanForce(6),MeanForce(4),. 

MeanForce(5)); 

fprintf(manyrowsfid,’\n'); 


Appendix  F.  Output  File  from  AnalyzemodXlsSF.m  called 
by  AutoanalyzeXls.m  for  Steady  Force  Tests 

This  is  part  of  the  output  file  from  AnalyzemodXlsSF.m  when  called  by 
AutoanalyzeXls.m  to  analyze  steady  force  tests.  The  two  rows  of  text  at  the  top  were 
added  by  the  author  for  clarification  and  are  not  actually  part  of  the  data  file.  When  using 
this  table  please  note  the  following: 

SF67-72,  80-90,  144-185  were  performed  with  control  surfaces  removed. 

SF73  was  a  test  to  determine  the  effect  of  system  noise  with  no  motion. 

SF74  was  performed  with  no  motion  but  with  a  4  oz  weight  attached  near  midships. 

SF75  was  performed  with  no  motion  but  with  a  2  oz  weight  attached  near  midships. 


Depth 

Subm. 

Spd 

Xdist 

Time 

I  Pitch  | 

HbuBI 

X 

Y 

Z 

Roll 

■KjcjEH 

■ 

■■ 

cm 

mm 

cm/s 

cm 

s 

|  Degrees 

N 

N 

N 

N-m 

SF126 

04/07/03 

79  5 

501.5 

100 

700 

10 

4 

0 

0 

0 

1.5582 

-0.3684 

-0.1921 

0.0073 

0.0015 

-0.0231 

SF127 

04/07/03 

79.5 

501.5 

75 

700 

13 

4 

0 

0 

0 

1.7731 

-0.0883 

-0.192 

0.0027 

0.0111 

-0.0039 

SF128 

04/07/03 

79.5 

356.2 

100 

700 

10 

8 

0 

0 

0 

-4.1236 

-0.1895 

-0.6621 

0.0044 

0.0651 

-0.0037 

SF129 

04/07/03 

79  5 

356.2 

75 

700 

13 

8 

0 

0 

0 

-2.5511 

-0.2603 

-0.4944 

0.0051 

0.0613 

-0  0105 

SF130 

04/07/03 

79.5 

210.9 

100 

700 

10 

8 

0 

0 

0 

5.4903 

0.0547 

-0.8964 

0.0277 

0.0674 

0.0062 

SF131 

04/07/03 

79.5 

210.9 

75 

700 

13 

8 

0 

0 

0 

8.7811 

0.239 

-0.7814 

0.026 

0.0676 

0,0145 

SF132 

04/07/03 

79.5 

501.5 

100 

700 

10 

8 

0 

0 

0 

-1.9805 

-0.1559 

-0.6271 

-0.0014 

0,0606 

-0.0079 

SF133 

04/07/03 

79.5 

501.5 

75 

700 

13 

8 

0 

0 

0 

1.3509 

-0.0493 

-0.4612 

0.0006 

0.0524 

-0.0012 

SF134 

04/07/03 

79  5 

356.2 

100 

700 

10 

12 

0 

0 

0 

2.9365 

-0.4546 

-1.0079 

0.0183 

0.1172 

-0.0088 

SF135 

04/07/03 

79.5 

356.2 

75 

700 

13 

12 

0 

0 

0 

0.1206 

-0.2504 

-0.7061 

0.0112 

0.0867 

-0.0135 

SF136 

04/07/03 

79.5 

210.9 

100 

700 

10 

12 

0 

0 

0 

10.111 

0.2624 

-1.6039 

0.027 

0.1288 

0.0077 

SF137 

04/07/03 

79.5 

210.9 

75 

700 

13 

12 

0 

0 

0 

13.6603 

0.2164 

-1.1318 

0.0324 

0.1141 

0.0176 

SF138 

04/07/03 

79.5 

501.5 

100 

700 

10 

12 

0 

0 

0 

1.0817 

-0  526 

-0.8771 

0.0042 

0.1005 

-0.0276 

SF139 

04/07/03 

79.5 

501.5 

75 

700 

13 

12 

0 

0 

0 

0.7022 

-0.1462 

-0.7014 

0.0016 

0.0688 

-0.004 

SF140 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

8 

0 

0 

8.2723 

08434 

0.7545 

-0.0361 

-0.0988 

0.194 

SF141 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

8 

0 

0 

6.0705 

0.7525 

0.0117 

-0.006 

-0.0367 

0.134 

SF142 

04/07/03 

79.5 

210.9 

100 

700 

10 

8 

0 

0 

0 

14.5032 

0.172 

-0.7051 

0.0024 

0.0661 

0.0052 

SF143 

04/07/03 

79.5 

210.9 

75 

700 

13 

4 

0 

0 

0 

19.6126 

0.2523 

-0.4773 

0.005 

0.013 

0.0121 

SF144 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

0 

0 

0 

11.9977 

-0.3885 

0.1854 

0.0021 

-0.0888 

-0.0355 

SF145 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

0 

0 

0 

12.6364 

0.0928 

-0.0909 

-0.0048 

-0.0267 

-0.0034 

SF146 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

0 

0 

0 

11.2486 

0.0564 

-0.3442 

0.0016 

-0.0728 

-0.0127 

SF147 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

0 

0 

0 

12.6827 

0.0567 

-0.2317 

0.004 

-0.0304 

-0.008 

SF148 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

0 

0 

0 

7.8708 

-0.3903 

0.2886 

-0.0143 

-0.0798 

-0.0368 

SF149 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

0 

0 

0 

7.3999 

-0.2902 

0.1767 

-0.0129 

-0.0526 

-0.0216 

SF150 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

4 

0 

0 

-0.8426 

0.0721 

0.4963 

-0.0322 

-0.0894 

0.1153 

SF151 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

4 

0 

0 

2.0863 

-0.1088 

0.6661 

-0.0248 

-0.0585 

0.0663 

SF152 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

4 

0 

0 

8.6837 

0.2023 

0.3061 

-0.0066 

-0.082 

0.1125 

SF153 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

4 

0 

0 

11.8808 

0.2309 

-0.4769 

-0.0046 

-0.003 

0.0779 

SF154 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

4 

0 

0 

4.0575 

0.1175 

0.3141 

-0.0248 

-0.0839 

0.1181 

SF155 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

4 

0 

0 

11.2072 

0.0949 

0.188 

-0.0183 

-0.0416 

0.0777 

SF156 

04/07/03 

795 

356.2 

100 

700 

10 

0 

8 

0 

0 

1.8573 

0.3133 

0.6259 

-0.0283 

-0.0693 

0.2394 

SF157 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

8 

0 

0 

10.5516 

0.6441 

-0.2885 

0.0098 

0.0026 

0.1728 

SF158 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

8 

0 

0 

12.9229 

1.0197 

0.0742 

-0.0048 

-0.0546 

0.2681 

SF159 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

8 

0 

0 

11.4304 

0.4812 

0.8981 

-0.0097 

-0.0339 

0.1663 

SF160 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

8 

0 

0 

-2.0303 

0.4854 

0.2242 

-0.0035 

-0.0645 

0.253 

SF161 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

8 

0 

0 

4.6134 

0.7048 

0.1588 

-0.0186 

-0.0291 

0.1856 

SF162 

04/07/03 

79.5 

356.2 

100 

700 

10 

0 

12 

0 

0 

-5.1568 

1.2942 

0.4418 

-0.0165 

-0.0632 

0.4016 

SF163 

04/07/03 

79.5 

356.2 

75 

700 

13 

0 

12 

0 

0 

3.8291 

0.9787 

-0.1571 

0.0085 

0.0067 

0.2613 

SF164 

04/07/03 

79.5 

210.9 

100 

700 

10 

0 

12 

0 

0 

10.272 

1.5353 

0.1432 

-0.0012 

-0.038 

0.4029 

SF165 

04/07/03 

79.5 

210.9 

75 

700 

13 

0 

12 

0 

0 

13.0686 

1.065 

-0.8274 

0.0037 

-0.0076 

0.2577 

SF166 

04/07/03 

79.5 

501.5 

100 

700 

10 

0 

12 

0 

0 

-7.5024 

1.4947 

0.18 

-0.017 

-0.0564 

0.4133 

SF167 

04/07/03 

79.5 

501.5 

75 

700 

13 

0 

12 

0 

0 

6.2886 

0.9665 

-0.0814 

0.0012 

-0.0278 

0.2643 

SF168 

04/07/03 

79.5 

356.2 

100 

700 

10 

4 

0 

0 

0 

7.7978 

-0.6456 

0.5741 

-0.0051 

0.0123 

-0.0521 

SF169 

04/07/03 

79.5 

356.2 

75 

700 

13 

4 

0 

0 

0 

9.444 

-0.1815 

0.1395 

-0.0051 

0.0376 

-0.0216 

SF170 

04/07/03 

79.5 

210.9 

100 

700 

10 

4 

0 

0 

0 

3.7602 

0.8067 

-0.3934 

-0.025 

0.0252 

0.0285 

SF171 

04/07/03 

79.5 

210.9 

75 

700 

13 

4 

0 

0 

0 

3.9603 

0.7468 

0.0696 

-0.0282 

0.0341 

0.04 

SF172 

04/07/03 

79.5 

501.5 

100 

700 

10 

4 

0 

0 

0 

-3.5378 

-0.1568 

-0.0781 

0.0035 

0.0379 

-0.027 

SF173 

04/07/03 

79.5 

501.5 

75 

700 

13 

4 

0 

0 

0 

-3.392 

-0.1565 

0.047 

0.0039 

00278 

-0.0167 

SF174 

04/07/03 

79.5 

356.2 

100 

700 

10 

8 

0 

0 

0 

8.2774 

0.2489 

-0.6666 

0.009 

0.1331 

-0.0013 

SF175 

04/07/03 

79.5 

356.2 

75 

700 

13 

8 

0 

0 

0 

15.9178 

0.3335 

-0.6135 

0.0198 

0.1163 

0.0029 

SF176 

•  04/07/03 

79.5 

210.9 

100 

700 

10 

8 

0 

0 

0 

14.9037 

0.3619 

-0.6767 

0.0128 

0.1374 

0.0167 

SF177 

04/07/03 

79.5 

210.9 

75 

700 

13 

8 

0 

0 

0 

19.9008 

0.7155 

-0.5469 

0.0154 

01086 

0.0298 

SF178 

04/07/03 

79.5 

501.5 

100 

700 

10 

8 

0 

0 

0 

13.905 

-0.3111 

-0.2381 

0.0055 

0.1145 

-0.0191 

SF179 

04/07/03 

79.5 

501.5 

75 

700 

13 

8 

0 

0 

0 

12.3352 

-0.1277 

-0.1732 

0.002 

0.0754 

-0.0116 

SF180 

04/07/03 

79.5 

3562 

100 

700 

10 

12 

0 

0 

0 

18.163 

0.2045 

-1.0493 

0.0147 

0.2235 

-0.003 

SF181 

04/07/03 

79.5 

356.2 

75 

700 

13 

12 

0 

0 

0 

17.4979 

0.1598 

-0.6601 

0.0119 

0.1675 

0.008 

SF182 

04/07/03 

79.5 

210.9 

100 

700 

10 

12 

0 

0 

0 

14.9671 

-0.0355 

-0.5625 

0.0014 

0.1805 

-0.0027 

SF183 

04/07/03 

79.5 

210.9 

75 

700 

13 

12 

0 

0  0 

16.0427 

-0.D323 

-0.2541 

0.003 

0.1056 

-0.0015 

SF184 

04/07/03 

795 

501.5 

100 

700 

10 

12 

0 

0 

0 

2.9632 

-0.1106 

-0.8265 

0.0007 

02131 

-0.0102 

SF185 

04/07/03 

795 

501.5 

75 

700 

13 

12 

0 

0 

0 

0.0831 

-0.2192 

-0.3446 

-0.005 

0.1249 

-0.0096 

# 


Appendix  G.  AnalyzemodXls.m 

%  AnalyzemodXls.m 
%  Erik  Oiler,  2003 

%  Series  Specific  Treatment 

%  Series  7:  Inertial  Results  are  not  subtracted  for  Mass  Matrix  Tests. 
%  Variables 

%  ActualPhasefl  ,f2,f3:  The  actual  phase  of  the  forces  at  the 
%  oscillation  frequency. 

%  Ampl:  The  amplitude  of  the  data  point. 

%  avg(13):  The  mean  value  of  the  force  and  motions  for  each 

%  experiment. 

%  B(6,6):  The  inverse  sensitivity  matrix  for  the  load  cell. 

%  C(6):  The  values  of  capacitances  in  the  electrical  filters  in 

%  microfarads. 

%  CF(6,6):  The  conversion  factor  matrix  for  the  load  cell. 

%  CFtemp(6):  The  conversion  factor  vector  for  the  load  cell  data. 
%  Data(NumPoints,15):  The  Datain  matrix  after  interpolation  for 
%  constant  time  intervals. 


%  Date: 

The  date  the  experiment  was  performed. 

%  DateNum: 

The  numeric  MATLAB  code  for  the  date  and  time  the 

%  experiment  was 

% 

performed. 

%  Datain: 

Holds  all  of  the  input  data  as  read  from  the  xls  file. 

% 

Column  1: 

Y  force 

% 

2 

Z  Force 

% 

3 

X  Force 

% 

4 

Pitch  Moment 

% 

5 

Yaw  Moment 

% 

6 

Roll  Moment 

% 

7 

Reserved  for  wave  data 

% 

8 

Reserved  for  wave  data 

% 

9 

X  position 

% 

10:  Y  position 

% 

1 

1 :  Z  Position 

% 

12:  Yaw  Position 

% 

1. 

3:  Pitch  Position 

% 

14:  Time 

% 

i 

5:  Date 

% 

16:  Time  Interval  between  samples 

%  DeltaT: 

The  average  time  interval  between  samples  in  seconds. 

%  Depth: 

The  depth  of  the  experiment  in  cm. 

%  dropgap: 

The  number  of  data  samples  dropped  from  the  start  of 

%  the  data  in  order  to  allow  for  acceleration. 

%  DrivingFreq:  The  frequency  of  driving  force  oscillation  in  Hz. 

%  This  is  omega  converted  into  Hz. 

%  Duration:  The  total  length  of  time  of  the  experiment. 

%  existencecheck:  Holds  the  results  of  the  check  for  the  input  file. 

%  fl :  The  ordered  frequency  of  oscillation  in  Hz. 

%  f2:  The  first  harmonic  of  the  oscillation  frequency  in  Hz. 

%  f3:  The  second  harmonic  of  the  oscillation  frequency  in  Hz. 

%  fname:  The  name  of  the  file  to  be  analyzed. 

%  Fpitch:  Calculated  inertial  force  in  pitch  in  Newton -meters. 
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%  Fx:  Calculated  inertial  force  in  the  X  direction  in  Newtons. 

%  Fy:  Calculated  inertial  force  in  the  Y  direction  in  Newtons. 

%  Fyaw:  Calculated  inertial  force  in  yaw  in  Newton-meters. 

%  Fz:  Calculated  inertial  force  in  the  Y  direction  in  Newtons. 

%  Gain:  The  gain  of  the  load  cell  amplifier. 

%  infname:  The  name  of  the  input  file  with  the  .xls  extension  added. 

%  indexhigh:  The  index  of  the  row  of  data  representing  the  frequency 

%  to  begin  filtering. 

%  InertialForce:  The  Calculated  Inertial  Force  in  N  or  N-m. 

%  Interval:  The  time  duration  of  the  fft  data. 

%  lyy:  First  moment  of  inertia  abouthte  Y-axis. 

%  Izz:  First  moment  of  inertia  about  the  Z-axis. 

%  k:  An  index  used  in  filtering. 

%  LocalMax:  The  local  maximum  of  data.  Used  in  finding  the  peak 

%  frequencies. 

%  LocalMaxCounter:  An  index  used  to  identify  the  lines  of  data  for  the 

%  local  maxima. 

%  look:  The  size  of  the  gap  between  frequency  intervals. 

%  manyrowsfid:  The  id  number  of  the  series  row  output  file. 

%  From  the  calling  file. 

%  mass:  The  mass  of  the  model  full  of  water  in  kg. 

%  MotionPhi:  The  phase  of  the  driving  force  motion  at  the  frequency 

%  of  the  local  maxima. 

%  nonzeroplanes:  From  the  calling  file.  Indicates  nonzero  stem 
%  planes  or  rudder. 

%  NumLines:  The  number  of  lines  in  the  input  file  containing  force 
%  and  motion  data. 

%  NumLines2:  Used  to  determine  when  the  model  stopped  traveling  in 
%  the  X  direction. 

%  NumPeriods:  The  number  of  periods  that  occur  over  the  duration  of 
%  the  experiment. 

%  NumPoints:  The  number  of  points  that  will  be  used  for  the  FFT. 

%  omega:  The  frequency  of  oscillation  in  radians  per  second. 

%  Period:  The  time  in  seconds  of  one  cycle  of  motion. 

%  Phaselndexfl  ,f2,f3:  The  line  of  Data  containing  the  actual  phase  of 

%  the  motion  at  the  oscillation  frequency. 

%  phi:  The  phase  shift  caused  by  dropping  the  first  1 .2  seconds  of 

%  data. 

%  PitchAmp:  Amplitude  of  Pitch  oscillation  in  degrees. 

%  PitchAngle:  Steady  pitch  angle.  Positive  is  nose  up. 

%  PitchDist:  Ordered  distance  of  travel  in  the  Pitch  direction  in 

%  degrees. 

%  pitchfactor:  The  control  factor  for  gantry  motion  in  pitch. 

%  PitchFreq:  Frequency  of  Pitch  oscillation  in  Hz. 

%  PitchPhase:  Ordered  phase  of  Pitch  oscillation  in  degrees. 

%  PitchRadius:  The  distance  from  the  center  of  gravity  to  the  axis  of 
%  pitch  in  cm. 

%  PitchSpd:  Speed  in  the  Pitch  direction  in  degreees/sec. 

%  PitchSupSpd:  Superimposed  speed  in  the  Pitch  direction  in  deg/sec. 

%  PitchSupDir:  Direction  of  superimposed  motion. 

%  R(6):  The  values  of  resistances  in  the  electrical  filters  in 

%  ohms. 

%  Results(3,24):  Holds  all  of  the  peaks  of  the  FFT. 

%  rowfname:  The  name  of  the  row  output  file. 

%  rowed:  The  file  id  number  of  the  row  output  file. 

%  Rudder:  The  angle  of  the  rudder  in  degrees. 


%  series:  The  series  of  data  to  be  analyzed. 

%  Comes  from  the  calling  program. 

%  1 :  Horizontal  Plane 

%  2:  Vertical  Plane 

%  3:  Pure  Sway 

%  4:  Pure  Heave 

%  5  :  Pure  Pitch 

%  6:  Pure  Yaw 

%  7:  Mass  Matrix 

%  8:  Steady  Force 

%  9:  Check  of  Inertial  Calculations 

%  10:  Miscellaneous 

%  11:  Oscillation  Tests 

%  SampleFreq:  Calculated  sample  frequency  based  on  average  interval 
%  between  datapoints. 

%  shiftlength:  The  distance  between  the  center  of  gravity  and  midships 
%  in  m. 

%  StemPlanes:  The  angle  of  the  stem  planes  in  degrees. 

%  Submergence:  The  submergence  of  the  model  to  the  top  of  the  hull  in 
%  mm.  From  the  calling  file. 

%  Summary (3, 24):  Holds  data  from  the  FFT  peak  closest  to  the  driving 
%  frequency. 

%  Time:  Ordered  length  of  time  of  the  experiment  in  seconds. 

%  TimeStart:  The  first  recorded  time  of  the  experiment. 

%  Vexc:  The  excitation  voltage  of  the  amplifier. 

%  write:  A  switch  for  whether  or  not  to  write  the  individual  row 
%  output  file. 


%  XAmp: 

%  XDist: 

%  xfactor: 

%  XFreq: 

%  Xg: 

% 

%  xhigh: 

%  XPhase: 

%  XSpd: 

%  XSupSpd: 

%  XSupDir: 

%  YAmp: 

%  Yaw  Amp: 

%  Yaw  Angle: 

%  stbd. 

%  YawDist: 

%  degrees. 


Amplitude  of  X  oscillation  in  cm. 

Ordered  distance  of  travel  in  the  X  direction  in  cm. 

The  control  factor  for  gantry  motion  in  the  X  direction. 
Frequency  of  X  oscillation  in  Hz. 

The  distance  from  the  model  center  of  gravity  to  the 
zero  point  of  the  load  cell. 

The  cutoff  frequency  for  low  pass  filtering. 

Ordered  phase  of  X  oscillation  in  degrees. 

Speed  in  the  X  direction  in  cm/sec. 

Superimposed  speed  in  the  X  direction  in  cm/sec. 
Direction  of  superimposed  motion. 

Amplitude  of  Y  oscillation  in  cm. 

Amplitude  of  Yaw  oscillation  in  degrees. 

The  steady  yaw  angle  in  degrees.  Positive  is  nose  to 


Ordered  distance  of  travel  in  the  Y aw  direction  in 


%  yawfactor: 

%  YawFreq: 

%  YawPhase: 
%  YawRadius: 
%  cm. 

%  YawSpd: 

%  YawSupSpd: 
%  YawSupDir: 
%  yfactor: 

%  YFreq: 

%  YPhase: 

%  YSpd: 


The  control  factor  for  gantry  motion  in  yaw. 

Frequency  of  Yaw  oscillation  in  Hz. 

Ordered  phase  of  Yaw  oscillation  in  degrees. 

The  distance  from  the  center  of  gravity  to  the  strut  in 

Speed  in  the  Yaw  direction  in  degreees/sec. 
Superimposed  speed  in  the  Yaw  direction  in  deg/sec. 
Direction  of  superimposed  motion. 

The  control  factor  for  gantry  motion  in  the  Y  direction. 
Frequency  of  Y  oscillation  in  Hz. 

Ordered  phase  of  Y  oscillation  in  degrees. 

Ordered  speed  of  travel  in  the  Y  direction  in  cm/sec. 


%  YSupSpd: 
%  YSupDir: 
%  ZAmp: 

%  ZDist: 

%  zfactor: 

%  ZFreq: 

%  ZPhase: 

%  ZSpd: 

%  ZSupSpd: 
%  ZSupDir: 


Superimposed  speed  in  the  Y  direction  in  cm/sec. 
Direction  of  superimposed  motion. 

Amplitude  of  Z  oscillation  in  mm. 

Ordered  distance  of  travel  in  the  Z  direction  in  mm. 
The  control  factor  for  gantry  motion  in  the  Z  direction. 
Frequency  of  Z  oscillation  in  Hz. 

Ordered  phase  of  Z  oscillation  in  degrees. 

Ordered  speed  of  travel  in  the  Z  direction  in  mm/sec. 

Superimposed  speed  in  the  Z  direction  in  mm/sec. 
Direction  of  superimposed  motion. 


%  Initialize  certain  variables. 

NumPoints  =  2048; 

Results  =  zeros(3,24);  %  Will  hold  all  of  the  Peaks  of  the  FFT. 

Summary  =  zeros(3,24);  %  Will  hold  the  FFT  peak  closest  to  the  driving 
%  frequency. 

[Data]  =  0; 

[Datain]=0; 


%  Get  the  input  data  and  find  the  length  of  the  file, 
infname  =  strcat(fhame,\xlsr); 

existencecheck  =  exist(infname);  %  Determines  if  the  input  file  exists, 
if  existencecheck  =0 

fprintf(’%s  does  not  existAn', infname) 
return 
end 

Datain  =  xlsread(infname);  %  Reads  the  input  file. 

NumLines  =  size(Datain,l)-4;  %  Gets  the  number  of  data  samples  in  the 
%  input  file. 


%  Get  the  date  from  the  input  file  and  convert  from  EXCEL  to 
%  MATLAB  format. 

Date=datestr(Datain(l,15)-36525,2); 

%  Build  a  date-time  string  to  ensure  the  correct  gains  are  applied. 
%  Add  693960  to  the  date  to  get  number  of  days  from  0000. 

%  Divide  the  number  of  minutes  by  86400  to  get  fractional  days. 
%  MATLAB  date  serial  numbers  are  in  the  form: 

%  days  since  0000. fraction  of  a  day 

DateN  um=Datain(  1 , 1 5)+693 960+Datain(  1 , 1 4)/86400; 


%  Record  input  motion  parameters. 
XAmp  =  Datain(NumLines+2,l); 
if  isnan(XAmp)=0 

XFreq  =  Datain(NumLines+2,2); 
XPhase  =  Datain  (NumLines+2, 3); 
XSupSpd  =  Datain(NumLines+2,4); 
XSupDir  =  Datain(NumLines+2,5); 
else 

XAmp=0; 

XFreq  -  0; 

XPhase  -  0; 

XSupSpd  =  0; 

XSupDir  =  0; 
end 


YAmp  =  Datain(NumLines+2?6); 


if  isnan(YAmp)=0 

YFreq  =  Datain(NumLines+2,7); 
YPhase  =  Datain(NumLines+2,8); 
YSupSpd  =  Datain(NumLines+2,9); 
YSupDir  =  Datain(NumLines+2,10); 
else 

YAmp=0; 

YFreq  =  0; 

YPhase  =  0; 

YSupSpd  =  0; 

YSupDir  =  0; 
end 

ZAmp  =  Datain(NumLines+2,l  1); 
if  isnan(ZAmp)=0 

ZFreq  =  Datain(NumLines+2,12); 
ZPhase  =  Datain(NumLines+2,13); 
ZSupSpd  =  Datain(NumLines+2J14); 
ZSupDir  =  Datain(NumLines+2,15); 
else 

ZAmp=0; 

ZFreq  =  0; 

ZPhase  =  0; 

ZSupSpd  =  0; 

ZSupDir  =  0; 
end 

YawAmp  =  Datain(NumLines+2, 1 6); 
if  isnan(YawAmp)=0 

YawFreq  =  Datain(NumLines+3,l); 
YawPhase  =  Datain(NumLines+3,2); 
YawSupSpd  =  Datain(NumLines+3?3); 
YawSupDir  =  Datain(NumLines+3,4); 
else 

YawAmp=0; 

YawFreq  =  0; 

YawPhase  =  0; 

YawSupSpd  =  0; 

YawSupDir  =  0; 
end 

PitchAmp  =  Datain(NumLines+3,5); 
if  isnan(PitchAmp)=0 
PitchFreq  =  Datain(NumLines+3,6); 
PitchPhase  =  Datain(NumLines+3,7); 
PitchSupSpd  =  Datain(NumLines+3,8) 
PitchSupDir  =  Datain(NumLines+3,9); 
else 

PitchAmp=0; 

PitchFreq  =  0; 

PitchPhase  =  0; 

PitchSupSpd  =  0; 

PitchSupDir  =  0; 
end 


XSpd  =  DatainfNumLines+3,10); 


if  isnan(XSpd)=l 
XSpd  -  0; 

if  ((isnan(XSupSpd)==0)&(XSupDir==  1 )) 

XSpd  =  XSupSpd; 
end 
end 

XDist  =  Datain(NumLines+3,l  1); 
if  isnan(XDist)=l 
XDist  =  0; 
end 

YSpd  =  Datain(NumLines+3,12); 

YDist  =  Datain(NumLines+3,l  3); 

ZSpd  =  Datain(NumLines+3,14); 

ZDist  =  Datain(NumLines+3,l  5); 

YawSpd  =  Datain(NumLines+3,16); 

YawDist  =  Datain(NumLines+4,l); 

PitchSpd  =  Datain(NumLines+4,2); 

PitchDist  =  Datain(NumLines+4,3); 

Time  =  Datain(NumLines+4,5); 

%  Determine  actual  sample  frequency  based  upon  the  measured  interval 
%  between  data  points.  Necessary  because  the  control  software  did  not 
%  always  sample  at  the  ordered  sample  rate. 

SampleFreq  =  l/mean(Datain(l:NumLines,16)); 

%  Determine  the  oscillation  frequency. 

fl=0; 

f2=0; 

0=0; 

if  (XFreq  ~=  0)&(isnan(XFreq)==0) 
fl  -  XFreq; 

£2  =  2*  XFreq; 

0  =  3*  XFreq; 
end 

if  ((YFreq  —  0)&(isnan(YFreq)==0)) 
fl  =  YFreq; 

0  =  2*  YFreq; 

0  =  3*  YFreq; 
end 

if  (ZFreq  ~=  0)&(isnan(ZFreq)==0) 
fl  =  ZFreq; 

£2  =  2*  ZFreq; 

0  =  3*  ZFreq; 
end 

if  (PitchFreq  ~  0)&(isnan(PitchFreq)=0) 
fl  =PitchFreq; 
f2  =  2*  PitchFreq; 

0  =  3*  PitchFreq; 


end 


# 


# 


$ 


» 


if  (YawFreq  ~  0)&(isnan(YawFreq)=0) 
fl  =  YawFreq; 
f2  =  2*  YawFreq; 

0  =  3*  YawFreq; 
end 

omega=fl*2*pi;  %  The  frequency  of  oscillation  in  radians  per  second. 

%  Drop  the  first  1 .2  seconds  of  data  to  account  for  acceleration.  At 
%  the  sample  frequency  of  25  Hz,  drop  tho  first  30  points.  Shift  the 
%  other  points  up  in  the  array. 

dropgap  =  round(l  .2*SampleFreq); 
for  I  =  1  :NumLines-dropgap; 

Datain(I,l  :16)=Datain(I+dropgap,l  :16); 
end 

NumLines  =  NumLines  -  dropgap; 

%  Drop  data  recorded  after  the  model  stopped  moving  along  the  X-axis 

%  for  non-Inertial  tests . 

if(XSpd“0) 

NumLines2==NumLines; 
fori  =  NumLines-l:-l:10 

if  (Datain(I,9)==Datain(HT  ,9))  %  Looks  for  constant  X  position. 

NumLines2  =  I; 
end 
end 

NumLines  —  NumLines2; 
end 

%  Make  the  time  of  good  data  be  an  integer  number  of  wavelengths. 

Period  =  1/fl; 

Duration  =  Datain(NumLines,  1 4)-Datain(  1,14); 

NumPeriods  =  Duration/Period; 

NumLines  =  round(SampleFreq*Period*fIoor(NumPeriods)); 

%  PREPROCESSING  STAGE 

%  Preprocess  the  input  file  for  analysis.  This  includes: 

%  1)  Accounting  for  the  phase  shift  caused  by  dropping  the  first  1 .2 
%  seconds. 

%  2)  Determining  the  mean  force  for  each  column  and  removing  that  mean 
%  from  the  column  data  to  get  dynamic  response, 

%  3)  Converting  force  data  from  voltages  to  forces  and  moments, 

%  4)  Converting  position  data  to  centimeters  for  x  and  y,  millimeters 
%  for  z,  and  degrees  for  the  angles, 

%  5)  Converting  the  time  past  midnight  to  time  of  data  run, 

%  6)  Finding  the  Mean  DeltaT, 

%  7)  Interpolating  the  data  and  time  to  produce  even  time  intervals 
%  and  the  associated  data. 

%  1,2)  Account  for  the  phase  shift  and  determine  the  mean  of  the  forces 
%  columns. 

for  J=  1:6 
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avg(J)  =  mean(Datain(l  :NumLines,J)); 
for  1  =  1  :NumLines 

Datain(l,J)  -  Datain(l,J)  -  avg(J); 
end 
end 

%  3,4)  Convert  force  data  from  voltages  to  forces  and  moments. 

%  Convert  position  data  to  centimeters  for  x  and  y,  millimeters  for 
%  z,  and  degrees  for  the  angles. 

TimeStart  =  Datain(l,14);  %  Used  for  converting  from  time  past  midnight 
%  to  time  of  run. 

%  The  inverse  sensitivity  matrix  (B) 

if  ((datenum(Date)>=datenum(2002,07, 1 7))&(datenum(Date)<datenum(2002, 1 0,3 1 ))) 
[B]  =  [  0.3901  0.0029 -0.0071  -0.0010 -0.0024 -0.0016; 

0.0023  0.3887  0.0016  0.0025-0.0039  0.0019; 

0.0139  0.0129  1.5006-0.0002-0.0215-0.0018; 

0.0000-0.0001  0.0016  0.0069  0.0000-0.0001; 

-0.0001  0.0000  0.0018  0.0000  0.0070  0.0000; 

-0.0003  -0.0002  0.0000  0.0000  0.0000  0.0108]; 
else  %valid  on  and  after  Oct  31,  2002 

[B]  =  [  0.3856  0.0020-0.0016-0.0008  0.0001  -0.0030; 

0.0023  0.3811  -0.0040  0.0012-0.0034  0.0031; 

0.0093  0.0024  1.5109  0.0068-0.0231  -0.0014; 

0.0001  0.0000  0.0013  0.0069  0.0000  0.0000; 

0.0000  0.0000  0.0007  0.0000  0.0069-0.0001; 

-0.0003  0.0001  0.0000  0.0000  0.0000  0.0106]; 

end 

%  Establish  Gains  and  Excitation  Voltage 
if  (DateNum>73 1434.61458)  %  August  6,  2002  1445. 

[Gain]  =  [4000;4000;1000;1000;1000;4000];  %  y  z  x  pitch  yaw  roll 
[Vexc]- [10;  10;  2.5;  10;  10;  10]; 
end 

%  Calculate  the  Conversion  Factors  (CF) 

CFtemp  =  Gain.*Vexc*10A-6; 

CF  =  zeros(6,6); 

CF(1,1)=  CFtemp(l); 

CF(2,2)=  CFtemp(2); 

CF(3,3)=  CFtemp(3); 

CF(4,4)=  CFtemp(4); 

CF(5,5)=  CFtemp(5); 

CF(6,6)=  CFtemp(6); 

%  Establish  multipliers 

%  yfactor  has  a  negative  multipler  to  convert  from  the  gantry  y  positive 

%  to  the  load  cell  y  positive  direction. 

if  ((DateNum>73 1 5 1 3.33333)&(DateNum<73 1551 .33333)) 

%  Between  Oct  24,  0800  and  Dec  1,  0800. 
xfactor  =  3850; 
yfactor  - -2410; 
zfactor  =  2114; 
pitchfactor  =  972; 


yawfactor=  1818; 

elseif  DateNum>73  3  55 1 .33333  %Dec  1,0800 
xfactor  =  3850; 
yfactor  =  -2410; 
zfactor  =  21 34; 
pitchfactor  =  155; 
yaw  factor  =  1818; 
else 

xfactor  =  3850; 
yfactor  = -2410; 
zfactor  =  4921; 
pitchfactor  =  1111; 
yawfactor=  1025; 
end 

%  Convert  voltages  to  forces  and  position  data  to  the  metric  system, 
for  1  =  1  :  NumLines 

Datain(I,l  :6)  =  (CFM  *B*Datain(I,l  :6)')'; 

Datain(I,9)  =  Datain(l,9)  /xfactor; 

Datain(I,10)  =  Datain(l,10)  /  yfactor; 

Datain(l,l  1)  =  Datain(I,l  1)  /  zfactor; 

Datain(I,12)  =  Datain(I,12)  /  yawfactor; 

Datain(I,13)  =  Datain(I,13)  /  pitchfactor; 

Datain(I,14)  =  Datain(l,14)  -  TimeStart; 
end 

%  Convert  the  position  data  to  position  of  midships  vice  position  of 
%  the  strut. 

distance  =  .1093;  %  The  distance  from  the  strut  to  midships. 

PitchRadius  =  .2576;  %  Length  of  pitch  arm 
for  1  =  1 :  NumLines 

%  Effect  of  yaw  on  y  position. 

Datain(I,10)  =  Datain(I,10)-distance*sin(Datain(I,12)*pi/180); 

%  Effect  of  pitch  on  x  position. 

Datain(l,9)  =  Datain(l,9)  +  PitchRadius*  sin(Datain(I,l  3)*pi/l  80); 

%  Effect  of  pitch  on  z  position. 

Datain(I,l  1)  =Datain(l,l  l)-PitchRadius*(l-cos(Datain(l,l  l)*pi/180)); 
end 

%  Account  for  the  transducer  being  forward  of  midships. 

%  Reference:  Marine  Hydrodynamics  by  J.N.  Newman,  Sect  6.2 
%  waveheight  =  Acos(kx-wt+arbitrary  phase  shift)  (eq  7  page  240) 

%  Let  x  =  0  at  the  midships  location. 

%  The  wave  is  travelling  in  the  positive  x  direction. 

%  The  transducer  is  at  a  negative  x  direction  from  midships. 

%  Assume  the  wave  is  travelling  in  the  positive  x-direction. 

%  Use  midships  as  the  zero  reference  location. 

%  Shift  the  origin  of  the  coordinate  system  from  the  origin  of  the  load 
%  cell  to  vessel  amidships, 
shiftlength  -  0.0522; 

wavenum  =  omegaA2/9.81 ; 
for  I  =l:NumLines 

Datain(l,4)  =  Datain(I,4)  -  Datain(I,2)  *  shiftlength;  %  pitch 
Datain(I,5)  =  Datain(l,5)  +  Datain(I,l)  *  shiftlength;  %  yaw 


%  6)  Determine  the  DeltaT  from  the  input  file. 

DeltaT  =  1/SampleFreq; 

%  7)  Interpolate  data  points  to  produce  even  time  intervals  and 
%  associated  data.  By  this  process  the  input  datain  matrix  will  be 
%  interpolated  into  the  data  matrix. 

%  Start  by  building  even  time  intervals  and  the  date  column, 
for  I  =  1  :NumPoints; 

%  Interpolate  as  if  the  sample  freq  was  25  hz. 

Data(I,14)  =  (1-1)*  .04; 

Data(1, 1 5)  =  Datain(  1,15); 
end 

%  Interpolate  each  column  in  turn, 
for  J  =  1:13 

Data(l  ,J)  =  Datain(l  ,J); 
for  I  =  2:NumLines; 
for  K  =  1  :NumLines; 

if  (Datain(K,  1 4)>Data(1 , 1 4)) 
break 
end 
end 

KS  =  K  -  1; 

if  (Data(1, 1 4)<=Datain(NumLines,  1 4)) 

Data(I,J)  -  Datain(KS,J)  +  (Data(I,14)-Datain(KS,14))*... 
(Datain(K,J)-Datain(KS,J))/(Datain(K,14)-Datain(KS,14)); 

else 

Data(I,J)  =  0.0; 
end 
end 
end 


%  Identify  the  driving  force, 
if  YFreq  ~  0  &  isnan(YFreq)=0 
DrivingForce  =  10; 
omega  =  YFreq*2*pi; 
elseif  ZFreq~=  0  &  isnan(ZFreq)==0 
DrivingForce  =11; 
omega  =  ZFreq*2*pi; 
elseif  YawFreq—  0  &  isnan(YawFreq)==0 
DrivingForce  =12; 
omega  =  YawFreq*2*pi; 
elseif  PitchFreq~0  &  isnan(PitchFreq)==0 
DrivingForce  =  13  ; 
omega  =  PitchFreq*2*pi; 
elseif  XFreq  ~  0  &  isnan(XFreq)==0 
DrivingForce  =  9; 
omega  =  XFreq*2*pi; 
else 

DrivingForce=9; 

omega=0; 


%  PROCESSING  STAGE 
for  J  =  10:13 

Data(l  :NumLines,jy=Data(l  :NumLines,J)-mean(Data(l  :NumLines,J)); 
end 


%  Round  all  data  up  to  2048  points.  Pad  the  excess  with  zero's. 
Data(NumLines+l  :NumPoints,l :  1 3)  =  0.0; 

%  Convert  the  time  column  to  a  frequency  column. 

Interval  =  .04  *  NumPoints; 

for  I  =  1 :  NumPoints 

Data(I,l  4)=(I-1  )/Interval; 
end 

%  For  each  column,  perform  Fast  Fourier  Transformation,  low  pass 
%  filtering,  and  identify  the  frequencies  and  amplitudes  of  the 
%  forces. 

%  Perform  the  Fast  Fourier  Transformation, 
for  J=  1:  13 

Data(:,J)  =  fft(Data(:,J)); 

%  Transform  the  FFT  Coefficient  to  account  for  padding, 
for  I  =  1  .NumPoints 

Data(I,J)  =  Data(l,J)*  NumPoints/NumLines; 
end 
end 


%  Change  the  phase  of  all  data  such  that  the  phase  is  relative  to 
%  the  driving  motion. 


for  I  =  1  :NumLines 
for  J=  1:13 

Data(I,J)  =  Data(l,J)*exp(-i*angle(Data(I,DrivingForce))); 
end 
end 

%  Account  for  ELECTRICAL  FILTERS. 

[R]  =  [145500  144100  147400  149200  191600  146600]; 

[C]  =  [.073  .1 14  .081  .099  .089  .075  ]*10A-6; 

for  J=  1 :6 

for  I  =  1  :NumLines 

Data(I,J)  =  Data(I,J)*exp(-i*atan(omega*C(J)*R(J))) ; 
end 
end 


%  Determine  the  actual  phase  of  the  motion  at  the  oscillation  frequency 
%  Find  the  index  of  the  frequency  of  oscillation. 
look=  .5*1/Interval; 

Phaselndexfl  =  find(((fl  dook)<=Data(:,14))&(Data(:,14)<(fl+look))); 
Phaselndexf2  =  find(((f2-look)<=Data(:,  1 4))&(Data(: ,  1 4)<(f2+look))); 


# 


PhaselndexO  =  fmd(((f3-look)<-Data(:7 1 4))&(Data( : ,  1 4)<(B+Iook))); 
for  J=9:13 

AcUialFreq(J)=Data(Phaselndexfl  ,14); 

%  ActualPhasefl  (J)=  angle(Data(Phaselndexfl ,  J))*  1 80/pi ; 

ActualPhasef2(J)-  angle(Data(Phaselndexf2,J))*l  80/pi; 
ActualPhasef3(J)=  angle(Data(PhaseIndexf3,J))*  1 80/pi; 
ActualMotion(J)=sqrt(2.*2.*((abs(Data(Phaselndexfl,J))/NumPoints)A2)); 
end 

^  omega  =  ActualFreq(DrivingForce)*2*pi; 


* 


%  Find  forces  at  the  frequencies  of  oscillation, 
for  J  =  1 :6 

Summaryfl(3*J-2)=:Data(PhaseIndexfl,14);  %  freq 

Summaryfl(3*J-l)=sqrt(2.*2.*((abs(Data(PhaseIndexfl  ,J))/NumPoints)A2));  %  ampl 
Summaryfl(3*J)=angle(Data(PhaseIndexfl  ,J))*  1 80/pi-ActualPhasefl  (DrivingForce);  %phase 

Summaryf2(3*J-2):=Data(PhaseIndexf2,14);  %  freq 

Summaryf2(3*J-l)=sqrt(2.*2.*((abs(Data(PhaseIndexf2,J))/NumPoints)A2));  %  ampl 
Summaryf2(3*J>=angle(Data(PhaseIndexf2,J))*  1 80/pi- ActualPhasef2(DrivingForce);  %phase 

Summaryf3(3  *  J-2)=Data(PhaseIndexf3 , 1 4);  %  freq 

Summaryf3(3*J-l)=sqrt(2.*2.*((abs(Data(Phaselndexf3,J))/NumPoints)A2));%  ampl 
Summaryf3(3*J)=angle(Data(PhaseIndexf3,J))*  1 80/pi-ActualPhasef3(DrivingForce);  %phase 
end 

%  Calculate  the  inertial  forces  in  mks  units  for  all  series  except 
%  mass  matrix  evaluation  (7). 
mass=  3.67; 

Xg=0.03; 

Gyradius  =  0.0764; 

YawRadius  =  0.07; 

if  series  ~  7  %  Calculates  inertial  force  for  all  except  mass  matrix 
%  (inertial)  tests. 

Iyy  =  0.092509; 

Izz  =  0.1006; 

XAmp  =  ActualMotion(9); 

YAmp  =  ActualMotion(lO); 

ZAmp  =  ActualMotion(l  1); 

PitchAmp  =  ActualMotion(13); 

Yaw  Amp  =  ActualMotion(12); 

Fx  =  -omegaA2*mass*(XAmp/100  +  PitchRadius*(PitchAmp*pi/180)); 

Fy  =  -omegaA2*(mass*(YAmp/100  +  ( Yaw  Amp  *pi/ 180)*  YawRadius)  +  mass*Xg*(YawAmp*pi/l  80)); 
Fz-  -omegaA2*(mass*(ZAmp/1000  +  (Pitch Amp *pi/ 180)  *  (YawRadius))- 
mass*Xg*(PitchAmp*pi/l  80)); 

Fpitch  =  -omegaA2*(-mass*Xg*(ZAmp/1000  +  (PitchAmp*pi/180)  *  (YawRadius-Xg))  +  lyy  * 
(PitchAmp*pi/l  80)); 

Fyaw  =  -omegaA2  *  (mass*Xg*(YAmp/100  -  (YawAmp*pi/180)*YawRadius)+Izz*(YawAmp*pi/l  80)); 
else 

Iyy  =  0; 
lzz=0; 

Fx=0; 

Fy=0; 

Fz=0; 

Fpitch  =  0; 
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%  Calculate  the  hydrodynamic  response  at  the  oscillation  frequency  and 
%  its  harmonics. 

DrivingFreq  =  omega/(2*pi); 


for  J=  1:6 
switch  J 
case  1 

InertialForee  =  Fy; 
case  2 

InertialForee  =  Fz; 
case  3 

InertialForee  =  Fx; 
case  4 

InertialForee  =  Fpitch; 
case  5 

InertialForee  =  Fyaw; 
case  6 

InertialForee  -  0; 
end 


Results(l,4*J-3)  =  Summary  fl  (3  *J~2);  %  oscillation  freq 
Results(l,4*J-2)  =  Summaryfl(3*J-l);  %  total  force 

Results(l,4*J-l)  =  imag(Summaryfl  (3*  J- 1  ))+real(  Summary  fl  (3  *  J- 1 ))  +  InertialForee;  %  force  without 
inertial  component 

Results(l,4*J)  =  Summaryfl(3*J);  %  phase 


Results(2,4*J-3)  =  Summaryf2(3*J-2);  %  first  harmonic  of  oscillation  freq 
Results(2,4*J-2)  =  Summaryf2(3*J-l);  %  total  force 
Results(2,4*J-l)  =  Summary£2(3*J-l);% 

Results(2,4*J)  =  Summaryf2(3*J);  %  phase 


Results(3,4*J-3)  -  Summaryf3(3*J-2);  %  second  harmonic  of  oscillation  freq 
Results(3,4*J-2)  =  Summaryf3(3*J-l);  %  total  force 
Results(3,4*J-l)  =  Summaryf3(3*J-l);  %  force  without  inertial  component 
Results(3,4*J)  =  Summaryf3(3*J);  %  phase 

end 

%  Condition  the  Results  array. 

%  Ensure  the  phase  angles  are  +-  180  degrees. 

%  Ensure  the  amplitudes  are  all  positive, 
for  I  =1:3 
for  J  =  1 :6 

if  Results(I,4*J-l)  <0 

Results(I,4*J-l)  =  abs(Results(I,4*J-l)); 

Results(I,4*J)  =  Results(I,4*J)  +  180; 
end 

if  Results(I,4*J)  >  180 
Results(I,4*J)  =  Results (1,4*1)  -  360; 
elseif  ResuIts(I,4*J)  <-180 

Results(I,4*J)  =  Results(I,4*J)  +  360; 


%  Print  to  a  common  row  output  file. 
fprintf(rnanyrowsfid,'%s  \t  %s\t  %6.2f  ,fname,Date,Depth); 

fyrintf(manyrowsfid  ^t%6  Af  \t%5 . 1  f  \t  %5 .2f\t%5 . 1  f ,  Submergence,  XSpd,XDist,  Duration); 
fprintf(manyrowsfid,’  \t  %2.0f  \t%2.0f\t  %3.1f\t  %3.  If, Pitch  Angle,  Y aw  AngIe,SternPlanes, Rudder); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f  \t  %4.1f  \t  %4.1f  \t%8.5f  \t  %7.4f  \t  %7.4f  \t 
%4. 1  f ,  ActualFreq(9),ActualMotion(9),XPhase,  ActualPhasefl  (9),Results(  1 ,9: 1 2)); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f\t  %7.4f  \t  %4.1f,ResuIts(2,  9:12)); 
fprintf(manyrowsfid,'  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(3,  9:12)); 
fprintf(manyrowsfid/  \t  %8.5f  \t  %7.4f  \t  %4.1f  \t  %4.1f  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t 
%4. 1  f  ,ActualFreq(  1 0),  ActualMotion(  1 0),YPhase,  ActualPhasefl  ( 1 0),Results(  1 , 1 :4)); 
fprint^manyrowsfid,'  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(2, 1:4)); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(3, 1 :4)); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f  \t  %4. If  \t  %4.1f  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t 
%4.1f,ActualFreq(l  l),ActualMotion(l  1),  ZPhase,  ActualPhasefl  (1  l),Results(  1,5: 8)); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(2,  5:8)); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f \t  %7.4f  \t  %4.1f,Results(3,  5:8)); 
fprintf(manyrowsfid,’  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f  ,Results(l,21 :24)); 
fprintf(manyrowsfid,’  \t  %8.5f \t  %7.4f  \t  %7.4f \t  %4.1f,Results(2,  21 :24)); 
fprintf(manyrowsfid,’  \t  %8.5f \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(3,  21 :24)); 
fprintfCmanyrowsfid,’  \t  %8.5f  \t  %7.4f\t  %4.1f  \t  %4.1f  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t 
%4. 1  f ,  ActualFreq(  1 3),ActualMotion(l  3),PitchPhase,  ActualPhasefl  ( 1 3),Results(  1,13:16)); 
fprintf(manyrowsfid/  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(2,  13:16)); 
fprintf(manyrowsfid,f  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f,Results(3,  13:16)); 
fprintf(manyrowsfid/  \t  %8.5f  \t  %7.4f  \t%4.1f\t  %4.1f  \t  %8.5f\t  %7.4f\t  %7.4f\t 
%4. 1  f , ActualFreq(  1 2),  ActualMotion(  1 2),  Y  awPhase,  ActualPhasefl  ( 1 2),Results(  1,17 :20)); 
fprintf(manyrowsfid,’  \t  %8.5f\t  %7.4f  \t  %7.4f\t  %4.1f,Results(2,  17:20)); 
fprintf(manyrowsfid,'  \t  %8.5f  \t  %7.4f  \t  %7.4f  \t  %4.1f  \n',Results(3,  17:20)); 


Appendix  H.  Results  of  Inertial  Force  Calculation  Checks 


Appendix  I.  Selected  Portion  of  the  Output  File  from  AnalyzemodXls.m  for  the  Analysis  of 
the  Oscillation  Test  Series 

The  first  two  rows  have  been  added  by  the  author  for  clarity. _ 
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Appendix  J.  CoeffSolver.m 

%CoeffSolver.m 


%  Erik  Oiler 
%  Spring  2003 

%  Solves  for  coefficients  using  linear  regression. 

%  Added  mass  and  damping  coefficients  are  of  the  form 
%  Coeff  =  aO  +  al  *L/subm  +  a3*Fr 
%  Restoring  forces  equations  are  of  the  form 
%  Coefficient=aO  +  al*L/subm  +  a2*L/submA2  +a3*Fr 
clear  all; 
close  all; 


%  Solve  for  Restoring  Force  Coefficients 

%  Solve  for  Yaw  Induced  Restoring  Force  (Yuv) 

YuvMatrix  =  [-5.368865729 

-4.853454974 

-2.329293369 

-3.514754719 

-2.848702047 

-1.614922323]; 


SteadyY  awMatrix=[  1 
1  1.720624859 

1  1.260063003 

1  2.711810649 

1  1.720624859 

1  1.260063003 


2.711810649 

2.960549904 

1.587758772 

7.353916996 

2.960549904 

1.587758772 


7.353916996  0.289533426 

0.289533426 

0.289533426 

0.386044568 

0.386044568 

0.386044568]; 


YuvCoeffMatrix=SteadyYawMatrix\YuvMatrix; 

PredY  uvMatrix=SteadyY  awMatrix*  Y  uvCoeffMatrix ; 
for  i  =  1  :size(PredYuvMatrix,l) 

YuvDiff(i)  =  abs((PredYuvMatrix(i)-YuvMatrix(i)))/YuvMatrix(i); 
end 

Y  uvrms=norm(Y  uvDiff)/sqrt(size(  Y  uvMatrix,  1 )); 


%  Solve  for  Pitch  Induced  Restoring  Force  (Zuw),  C53 

ZuwMatrix  =  [1 .023955274 

2.269136274 

1.141877357 

1.220249361 

2.612356757 

1.801033337]; 


SteadyPitchMatrix=[l 

2.711810649 

7.353916996  0.289533426 

1 

1.720624859 

2.960549904 

0.289533426 

1 

1.260063003 

1.587758772 

0.289533426 

1 

2.714285714 

7.367346939 

0.386044568 

1 

1.718592965 

2.953561779 

0.386044568 

1 

1.259668508 

1.586764751 

0.386044568]; 

ZuwCoeffMatrix=SteadyPitchMatrix\ZuwMatrix; 


PredZuwMatrix=SteadyPitchMatrix*ZuwCoeffMatrix; 
for  i  =  l:size(PredZuwMatrix,l) 

ZuwDiff(i)  =  abs((PredZuwMatrix(i)-ZuwMatrix(i)))/ZuwMatrix(i); 
6  end 

Zuwrms=norm(ZuwDiff)/sqrt(size(ZuwMatrix,  1 )); 


%  Solve  for  Yaw  Induced  Restoring  Moment  (Nuv) 


NuvMatrix  =  [-1.020609587 

-0.996711818 

-0.682640735 

-0.800283816 

-0.737734321 

-0.586442502]; 


SteadyY  awMomMatrix=[  1 
1  1.720624859 

1  1.260063003 

1  2.711810649 

1  1.720624859 

1  1.260063003 


2.711810649  7.353916996 


2.960549904 

1.587758772 

7.353916996 

2.960549904 

1.587758772 


0.289533426 

0.289533426 

0.386044568 

0.386044568 

0.386044568]; 


0.289533426 


9  NuvCoeffMatrix=SteadyYawMomMatrix\NuvMatrix; 

PredNuvMatrix=SteadyYawMomMatrix*NuvCoeffMatrix; 
for  i  =  1  :size(PredNuvMatrix,l) 

NuvDiff(i)  =  abs((PredNuvMatrix(i)-NuvMatrix(i)))/NuvMatrix(i); 
end 

Nuvrms=norm(NuvDif!)/sqrt(size(NuvMatrix,  1 )); 


%  Solve  for  Pitch  Induced  Restoring  Moment  (Muw) 
MuwMatrix  =  [-0.626004966 


-0.594465093 

-0.71149292 

-0.401198364 

-0.619156216 

-0.599156406]; 


SteadyPitchMomMatrix=[  1 


2.711810649  7.353916996 


1.720624859 

1.260063003 

2.714285714 

1.718592965 

1.259668508 


2.960549904 

1.587758772 

7.367346939 

2.953561779 

1.586764751 


0.289533426 

0.289533426 

0.386044568 

0.386044568 

0.386044568]; 


0.289533426 


MuwCoeffMatrix=SteadyPitchMomMatrix\MuwMatrix; 

PredMuwMatrix=SteadyPitchMomMatrix*MuwCoeffMatrix; 

®  for  i  =  1  :size(PredMuwMatrix,l) 

MuwDiff(i)  =  abs((PredMuwMatrix(i)-MuwMatrix(i)))/MuwMatrix(i); 
end 

Muwrms=norm(MuwDiff)/ sqrt(size(MuwM  atrix,  1 )) ; 


%  Solve  for  Sway  induced  added  masses 


%  Solve  for  Yvdot 
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YvdotMatrix=[“0.017470 

-0.018969 

-0.016666 

-0.019724 

-0.017526]; 


SwayMatrix=[  1  1 .277450258  0. 1 276785 1 2 

1  2.749405234  0.383418954 

1  1.744402516  0.255740442 

1  1.277450258  0.383418954 

1  2.749405234  0.127678512]; 

YvdotCoeffMatrix  =  SwayMatrix\YvdotMatrix; 
PredY  vdotMatrix=SwayMatrix*  Y  vdotCoeffMatrix; 


sumofsquares=0; 


for  i  =  1  :size(PredY  vdotMatrix ,  1 ) 

YvdotDiff(i)  =  (PredYvdotMatrix(i)-YvdotMatrix(i))/YvdotMatrix(i); 
end 

Yvdotrms=norm(YvdotDif!)/sqrt(size(YvdotMatrix,l)); 

%  Solve  for  Yv 

YvMatrix=[-0.058883 

-0.027836 

-0.062066 

-0.027142 

-0.058392]; 

YvCoeffMatrix  =  SwayMatrix\YvMatrix; 

PredY  vMatrix=SwayMatrix*YvCoeffMatrix; 

sumofsquares=0; 

for  i  =  1 :  size(Pred  Y  vMatrix,  1 ) 

YvDiff(i)  -  abs((PredYvMatrix(i)-YvMatrix(i)))ArvMatrix(i); 
end 

Y  vnns=norm(Y  vDiff)/sqrt(size(Y  vMatrix,  1 )); 


%  Solve  for  Nvdot 

NvdotMatrix=[-0.000809 

-0.000692 

-0.000677 

-0.000862 

-0.000833]; 

NvdotCoeffMatrix  -  SwayMatrix\N vdotMatrix; 
PredNvdotMatrix=SwayMatrix*NvdotCoeffMatrix; 

for  i  =  1  :size(PredNvdotMatrix,  1) 

NvdotDiff(i)=abs(PredNvdotMatrix(i)-NvdotMatrix(i))/NvdotMatrix(i); 

end 

Nvdotrms=norm(NvdotDiff)/sqrt(size(NvdotMatrix,l)); 

%  Solve  for  Nv 

NvMatrix=[-0.004323 

-0.008228 


-0.005161 

-0.008666 

-0.004136]; 

NvCoeffMatrix  =  SwayMatrixYNvMatrix; 
PredNvMatrix=SwayMatrix*NvCoeffMatrix; 

for  i  =  l:size(PredNvMatrix,l) 

NvDiff(i)  =  (PredNvMatrix(i)-NvMatrix(i))/NvMatrix(i); 

end 

Nvrms=norm(NvDiff)/sqrt(size(NvMatrix,l)); 

%  Solve  for  coefficients  in  heave 
%  Solve  for  Zwdot 

ZwdotMatrix^-0.0 19162 
-0.022988 
-0.016744 
-0.020764 
-0.021686]; 


HeaveMatrix=[l 

1.421 

1 

2.547 

0.383 

1 

1.744 

0.256 

1 

1.421 

0.383 

1 

2.547 

0.128]; 

ZwdotCoeffMatrix  =  HeaveMatrixYZwdotMatrix; 
PredZwdotMatrix=HeaveMatrix*ZwdotCoeffMatrix; 

for  i  =  1  :size(PredZwdotMatrix,l) 

ZwdotDiff(i)=abs((PredZwdotMatrix(i)-ZwdotMatrix(i)))/ZwdotMatrix(i); 

end 

Zwdotrms=norm(ZwdotDiff)/sqrt(size(ZwdotMatrix,  1 )); 

%  Solve  for  Zw 

ZwMatrix=[-0.062616 

-0.029098 

-0.065158 

-0.024766 

-0.063094]; 

ZwCoeffMatrix  =  HeaveMatrix\ZwMatrix; 
PredZwMatrix=HeaveMatrix*ZwCoeffMatrix; 

for  i  =  1  :size(PredZwMatrix,l) 

ZwDiff(i)  =  abs((PredZwMatrix(i)-ZwMatrix(i)))/ZwMatrix(i); 
end 

Zwrms=norm(ZwDiff)/ sqrt(size(ZwMatrix,  1 )); 


%  Solve  for  Mwdot 

MwdotMatrix=[0.002835 

0.002002 

0.002847 


0.000781 

0.003218]; 

MwdotCoeffMatrix  =  HeaveMatrix\MwdotMatrix; 
PredMwdotMatrix=HeaveMatrix*MwdotCoeffMatrix; 

for  i  =  1  :size(PredMwdotMatrix,l) 

MwdotDiff(i)  =  abs((PredMwdotMatrix(i)-MwdotMatrix(i)))/MwdotMatrix(i); 
end 

Mwdotrms=norm(MwdotDiff)/sqrt(size(MwdotMatrix,l)); 

%  Solve  for  Mw 

MwMatrix=[0.0 1 845 1 

0.008741 

0.009078 

0.010553 

0.012233]; 

MwCoeffMatrix  =  HeaveMatrix\MwMatrix; 
PredMwMatrix=HeaveMatrix*MwCoeffMatrix; 

fori  =  l:size(PredMwMatrix,l) 

MwDiff(i)  =  abs((PredMwMatrix(i)-MwMatrix(i)))/MwMatrix(i); 
end 

Mwrms=norm(MwDiff)/sqrt(size(MwMatrix,l)); 


%  Solve  for  coefficients  in  yaw 
%  Solve  for  Yrdot 


Y  rdotMatrix=[0 .00 1338 

0.002209 

0.003194 

0.002134 

0.001086]; 


YawMatrix-[l 
1  1.277 

1  1.744 

1  2.749 

1  1.277 


2.749  0.128 

0.128 

0.256 

0.383 

0.383]; 


YrdotCoeffMatrix  =  YawMatrix\YrdotMatrix; 
PredYrdotMatrix^YawMatrix*  YrdotCoeffMatrix; 

for  i  =  l:size(PredYrdotMatrix,l) 

YrdotDiff(i)=abs((PredYrdotMatrix(i)-YrdotMatrix(i)))ArrdotMatrix(i); 

end 

Yrdotrms=norm(YrdotDiff)/sqrt(size(YrdotMatrix,l)); 


%  Solve  for  Yr 
Y  rMatrix-[0.02 1030 
0.018478 
0.015644 
0.014240 


0.016104]; 


YrCoeffMatrix  =  YawMatrix\YrMatrix; 
Pred  Y  rM  atrix= Y  awM  atrix  *  Y  rCoeffM  atrix ; 


fori=  l:size(PredYrMatrix,l) 

YrDiff(i)  =  abs((PredYrMatrix(i)-YrMatrix(i)))Af  rMatrix(i); 
end 

Yrrms=norm(YrDiff)/sqrt(size(YrMatrix?l)); 

%  Solve  for  Nrdot 

NrdotMatrix=[-0.00087 1 

-0.000844 

-0.00113666 

-0.001051 

-0.000960]; 

NrdotCoeffMatrix  =  YawMatrix\NrdotMatrix; 
PredNrdotMatrix=YawMatrix*NrdotCoeffMatrix; 

for  i  =  1  :size(PredNrdotMatrix,l) 

NrdotDiff(i)=abs((PredNrdotMatrix(i)-NrdotMatrix(i)))/NrdotMatrix(i); 

end 

Nrdotrms=norm(NrdotDiff)/sqrt(size(NrdotMatrix,l)); 

%  Solve  for  Nr 

NrMatrix=[-0.008514 

-0.005804 

-0.003476 

-0.004136 

-0.003916]; 

NrCoeffMatrix  =  YawMatrix\NrM  atrix; 
PredNrMatrix=YawMatrix*NrCoeffMatrix; 

for  i  =  1  :size(PredNrMatrix,l) 

NrDiff(i)  =  abs((PredNrMatrix(i)-NrMatrix(i)))/NrMatrix(i); 
end 

Nrrms=norm(NrDiff)/sqrt(size(NrMatrix,l)); 

%  Solve  for  coefficients  in  Pitch 
%  Solve  for  Zqdot 

ZqdotMatrix=[-0.00 1 497 1 96 

-0.000740791 

0.0036744]; 

PitchMatrix=[l  2 .749405234  0. 1 276785 1 2 

1  1.277450258  0.127678512 

1  2.749405234  0.383418954]; 

ZqdotCoeffMatrix  =  PitchMatrix\ZqdotMatrix; 
PredZqdotMatrix=PitchMatrix*ZqdotCoeffMatrix; 
for  i  =  1  :size(PredZqdotMatrix,l) 

ZqdotDiff(i)=abs((PredZqdotMatrix(i)-ZqdotMatrix(i)))/ZqdotMatrix(i); 


end 

Zqdotrms=norm(ZqdotDiff)/sqrt(size(ZqdotMatrix?l)); 


£  %  Solve  for  Zq 

ZqMatrix=[-0.0435061 32 

-0.048655759 

-0.009306293]; 


ZqCoeffMatrix  -  PitchMatrix\ZqMatrix; 
PredZqMatrix=PitchMatrix*ZqCoeffMatrix; 

fori  =  l:size(PredZqMatrix,l) 

ZqDiff(i)  =  abs((PredZqMatrix(i)-ZqMatrix(i)))/ZqMatrix(i); 
end 

Zqrms=norm(ZqDiff)/sqrt(size(ZqMatrix,  1 )); 

%  Solve  for  Mqdot 
MqdotMatrix=[-0.00082561 8 
-0.00098228 
-0.003345578]; 

MqdotCoeffMatrix  =  PitchMatrixVMqdotMatrix; 
PredMqdotMatrix=PitchMatrix*MqdotCoeffMatrix; 

for  i  =  1  :size(PredMqdotMatrix,l ) 

MqdotDiff(i)=abs((PredMqdotMatrix(i)-MqdotMatrix(i)))/MqdotMatrix(i); 

end 

Mqdotrms=norm(MqdotDiff)/sqrt(size(MqdotMatrix,l)); 


%  Solve  for  Mq 

MqMatrix=[-0.0085 13669 

-0.007481107 

-0.000642713]; 

MqCoeffMatrix  =  PitchMatrixYMqMatrix; 

PredMqMatrix=PitchMatrix*MqCoeffMatrix; 

for  i  =  1  :size(PredMqMatrix,l) 

MqDiff(i)  =  abs((PredMqMatrix(i)-MqMatrix(i)))/MqMatrix(i); 
end 

Mqrms=norm(MqDiff)/sqrt(size(MqMatrix,  1 )); 

%  Build  Output  File 
warning  off 
deleteCCoeffSolver.txf) 
warning  on 

fileid  =  fopenCCoeffSolver.txtVa'); 

fprintf(fileidJ’%s\t%s\t%s\t%s\t%s\t%s\t\nVCoeff,,lVL/Subm?>,(L/Subm)A2,,,FrVrms?); 
fprintf(fileid/%s\t%8.6f\t%8.6f\t  %8.6f  \t%8.6f\t%g  W/Yuv’^uvCoeffMatrix^uvrms); 
fprintf(Fileid/%s  \t  %8.6f  \t  %8.6f\t  %8.6f\t  %8.6f\t%g  W/Zuw\ZuwCoeffMatrix,Zuwrms); 
fprintf(fileid/%s  \t  %8.6f  \t  %8.6f\t  %8.6f  \t  %8.6f  \t%g  Xn’/Nuv^NuvCoeffMatrix^uvrms); 
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fprintf(fileid,'%s  \t  %8.6f  \t  %8.6f\t  %8.6f  \t%8.6f  \t  %g  \n’,'Muw',MuwCoeffMatrix,Muwrms); 
fprintf(fileid,'%s  \t%s  \t%s  \t%s  \t%s  \t\n’,’Coeff ,’l  ,,TL/Subm',’Fr,,,rms’); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’,’Yrdot’,YrdotCoeffMatrix,Yrdotnns); 
fprintf(fileid,'%s  \t  %8.6f  \t  %8,6f\t  %8.6f\t  %g  \n','Yr’,YrCoeffMatrix,Yrrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’,'Y vdot’.Y vdotCoeffMatrix,Y vdotrms); 
fprintf(fileid,'%s  \t  %8.6f  \t  %8.6f\t  %8.6f\t  %g  \nYYv',YvCoeffMatrix,Yvrms); 
fprintf(fileid,'%s  \t  %8.6f  \t  %8.6f\t  %8.6f\t  %g  \n',’Zqdot’,ZqdotCoefMatrix,Zqdotrms); 
fprintf(fileid,'%s  \t  %8.6f  \t  %8.6f\t  %8.6f\t  %g  \nyZq',ZqCoeffMatrix,  Zqrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n',’Zwdot’,ZwdotCoeffMatrix,Zwdotrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’,’Zw',ZwCoeffMatrix,  Zwrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’,'Mqdot’,MqdotCoef£Matrix,Mqdotrms); 
fprintf(fileid,'%s  \t  %8.6f  \t  %8.6f\t  %8.6f\t  %g  \n’,’Mq',MqCoeffMatrix,  Mqrms); 
fprintf(fileid,’%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’,’Mwdot’,MwdotCoef{Matrix,Mwdotnns); 
fprintf(fileid,'%s  \t  %8.6f  \t  %8.6f\t  %8.6At  %g  \n’,'Mw’,MwCoefifMatrix,  Mwrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n','Nrdot',NrdotCoeffMatrix,  Nrdotrms); 
fprintf(fileid,’%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’;Nr’,NrCoeffMatrix,  Nrrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n’,'Nvdot’,NvdotCoeffMatrix,  Nvdotrms); 
fprintf(fileid,'%s  \t  %8.6f \t  %8.6f\t  %8.6f\t  %g  \n',’Nv’,NvCoeffMatrix,  Nvrms); 
status  =  fclose(fileid); 


Appendix  K.  Output  of  CoeffSolver.m 

This  file  shows  the  output  of  the  file  CoeffSolver.m.  To  derive  the  equation  for 
coefficient  A ' ,  find  the  sum  of  the  products  of  the  elements  in  the  row  labeled  A  and  the 
value  of  the  header  for  the  column  of  the  element,  excepting  the  last  column.  In 
mathematical  form,  the  equation  looks  is 

A'  =  y, (A, Header,)  (44) 

j= i 

where  n  =4  for  restoring  force  coefficients  and  n=3  for  all  other  coefficients.  The  root 
mean  square  of  the  difference  between  the  predicted  and  measured  values  are  shown  in 
the  final  column. 


Coeff 

1 

L/Subm 

(L/Subm)A2 

Fr 

rms 

Yuv 

3.035943 

-11.2325 

2.39970 

15.79519 

0.132812 

Zuw 

-7.723764 

9.156976 

-2.36537 

4.182974 

0.0752513 

Nuv 

0.020017 

-1.45297 

0.31796 

1.987686 

0.050452 

Muw 

-1.168727 

0.127681 

-0.00757 

1.079278 

0.0951314 

Coeff 

1 

L/Subm 

Fr 

rms 

Yrdot 

0.002324 

-8.7E-05 

-0.00063 

0.462378 

Yr 

0.020906 

0.000403 

-0.01798 

0.072234 

Yvdot 

-0.016345 

0.000061 

-0.00722 

0.043182 

Yv 

-0.081458 

0.001774 

0.12175 

0.148568 

Zqdot 

-0.002666 

-0.00051 

0.02022 

2.00E-16 

Zq 

-0.070199 

0.003498 

0.13373 

8.66E-16 

Zwdot 

-0.013633 

-0.00268 

-0.00567 

0.088254 

Zw 

-0.086685 

0.00091 

0.14075 

0.170412 

Mqdot 

0.00014 

0.000106 

-0.00985 

1.66E-16 

Mq 

-0.010515 

-0.0007 

0.03078 

1.18E-15 

Mwdot 

0.002825 

0.000594 

-0.00641 

0.271656 

Mw 

0.023224 

-0.00294 

-0.02236 

0.190629 

Nrdot 

-0.000786 

-1.9E-05 

-0.00058 

0.080731 

Nr 

-0.006013 

-0.00117 

0.01230 

0.232831 

Nvdot 

-0.00089 

0.000037 

0.00017 

0.092017 

Nv 

-0.00207 

0.000093 

-0.01649 

0.093355 

Appendix  L.  Model  Scale  Experiments  Performed  at  MIT  to 
Determine  the  Restoring  Forces  and  Moments  due  to  Body 
Angle 


Test  # 

Water 

Depth 

Subm. 

Yaw 

Pitch 

Velocity 

in 

m 

Degrees 

Degrees 

m/sec 

SF98 

0.795 

0.3962 

0 

0 

1 

SF99 

0.795 

0.3962 

0 

0 

0.75 

SF100 

0.795 

0.2509 

0 

0 

1 

SF101 

0.795 

0.2509 

0 

6 

0.75 

SF102 

0.795 

0.5415 

0 

0 

1 

SF103 

0.795 

0.5415 

0 

0 

0.75 

SF104 

0.795 

0.3962 

4 

0 

1 

SF105 

0.795 

0.3962 

4 

0 

0.75 

SF106 

0.795 

0.2509 

4 

0 

1 

SF107 

0.795 

0.2509 

4 

0 

0.75 

SF108 

0.795 

0.5415 

4 

0 

1 

SF109 

0.795 

0.5415 

4 

0 

0.75 

SF110 

0.795 

0.3962 

8 

0 

1 

SFlll 

0.795 

0.3962 

8 

0 

0.75 

SF112 

0.795 

0.2509 

8 

0 

1 

SF113 

0.795 

0.2509 

8 

0 

0.75 

SF114 

0.795 

0.5415 

8 

0 

1 

SF115 

0.795 

0.5415 

8 

0 

0.75 

SF116 

0.795 

0.3962 

12 

0 

1 

SF117 

0.795 

0.3962 

12 

0 

0.75 

SF118 

0.795 

0.2509 

12 

0 

1 

SF119 

0.795 

0.2509 

12 

0 

0.75 

SF120 

0.795 

0.5415 

12 

0 

1 

SF121 

0.795 

0.5415 

12 

0 

0.75 

SF122 

0.795 

0.3962 

0 

4 

1 

SF123 

0.795 

0.3962 

0 

4 

0.75 

SF124 

0.795 

0.2509 

0 

4 

1 

SF125 

0.795 

0.2509 

0 

4 

0.75 

SF126 

0.795 

0.5415 

0 

4 

1 

SF127 

0.795 

0.5415 

0 

4 

0.75 

SF128 

0.795 

0.3962 

0 

8 

1 

SF129 

0.795 

0.3962 

0 

8 

0.75 

SF130 

0.795 

0.2509 

0 

8 

1 

SF131 

0.795 

0.2509 

0 

8 

0.75 

SF132 

0.795 

0.5415 

0 

8 

1 

SF133 

0.795 

0.5415 

0 

8 

0.75 

SF134 

0.795 

0.3962 

0 

12 

1 

SF135 

0.795 

0.3962 

o" 

12 

0.75 

SF136 

0.795 

0.2509 

0 

12 

1 

SF137 

0.795 

0.2509 

0 

12 

0.75 

SF138 

0.795 

0.5415 

0 

12 

1 

SF139 

0.795 

0.5415 

0 

12 

0.75 

100 


/ 


Appendix  M.  Results  of  Model  Scale  Experiments  Performed 
at  MIT  to  Determine  the  Restoring  Forces  and  Moments  due  to 


Body  Angle 

Shaded  cells  represent  values  greater  than  1.15  standard  deviations  away  from  the  mean 


Test# 

Subm. 

Spd 

Pitch 

Yaw 

X 

Y 

z 

L_ _ 

Roll 

Pitch 

Yaw 

Normalized 
Sway  Force 

Normalized 
Heave  Force 

Normalized 
Yaw  Moment 

Normalized 
Pitch  Moment 

SF101 

252 

75 

0 

0 

EEBH 

EEBsH 

BfiligliH 

cheese 

-0.0368 

0 

6 

0 

SF107 

252 

75 

0 

4 

■he  raw 

Bifchikl 

ebhe^ 

0.0647 

heeded 

HEE5ED3 

SF113 

252 

75 

0 

8 

■HttdFI 

MtMkkl 

-0.033 

0.1305 

EEEH^i 

SF119 

252 

75 

0 

12 

HEED 

K M 

EE 

MOUKEil 

0.0022 

1.799 

^ES^ 

H  H 

SF125 

252 

75 

4 

0 

■«es*i 

0.0337 

0.0086 

keih 

EHEDZ£E] 

SF143 

252 

75 

4 

0 

HEDDD3 

0.013 

0.0121 

bheddes 

HHDSQ3 

SF131 

252 

75 

8 

0 

HI 

kekeeb 

0.0145 

ehesees 

SF137 

252 

75 

12 

0 

BE 

BHHE3 

0.1141 

0.0176 

^kbisb 

0.1557 

SF99 

398 

75 

0 

0 

Mdhikld 

-0.0474 

EBNHBE] 

0 

^H^Hj 

HIHHBu] 

SF105 

398 

75 

0 

4 

HESEa 

keep 

ITtEEEl 

BHEDDDS 

0.0572 

HEBSEQ 

hksezs 

SF111 

398 

75 

0 

8 

■i  fcbfeM 

ebbed 

-0.044 

^^^EHED] 

0.9341 

heedes 

SF117 

398 

75 

0 

12 

giOTa 

BTTi’dlH 

-0.0279 

HHEXESO 

1.5408 

BBEHE3 

SF123 

398 

75 

4 

0 

ebbed 

EEEE 

-0.003 

-0.0077 

-0.184 

flBBEDIEc] 

SF129 

398 

75 

8 

0 

EEESSSB 

eesd 

0.0613 

KHE££ES 

-0.719 

HiH 

SF135 

398 

75 

12 

0 

Belial 

■dh»H 

0.0867 

-0.0135 

-0.9307 

HH3EEE3 

SF103 

543 

75 

0 

0 

EESEB 

EH 

^WaMiliLn 
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Units  for  Interpreting  the  Results  of  Model  Scale  Experiments  Performed  at  MIT  to  Determine  the 
Restoring  Forces  and  Moments  due  to  Body  Angle 


Quantity 

Units 

Depth 

cm 

Submergence 

mm 

Speed 

Cm/s 

Angle 

Degrees 

Forces 

Newtons 

Moments 

N-m 
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